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We investigate worlds that lie on a slice through the parameter space of the Standard Model over 
which quark masses vary. We allow as many as three quarks to participate in nuclei, while fixing 
the mass of the electron and the average mass of the lightest baryon flavor multiplet. We classify 
as congenial worlds that satisfy the environmental constraint that the quark masses allow for stable 
nuclei with charge one, six, and eight, making organic chemistry possible. Whether a congenial world 
actually produces observers capable of measuring those quark masses depends on a multitude of 
historical contingencies, beginning with primordial nucleosynthesis and including other astrophysical 
processes, which we do not explore. Such constraints may be independently superimposed on our 
results. Environmental constraints such as the ones we study may be combined with information 
about the a priori distribution of quark masses over the landscape of possible universes to determine 
whether the measured values of the quark masses are determined environmentally, but our analysis 
is independent of such an anthropic approach. 

We estimate baryon masses as functions of quark masses via first-order perturbation theory in 
flavor SU{?>) breaking. We estimate nuclear masses as functions of the baryon masses using two 
separate tools: for a nucleus made of two baryon species, when possible we consider its analog in 
our world: a nucleus with a similar binding energy, up to Coulomb contributions. For heavy nuclei 
or nuclei made of more than two baryons, we develop a generalized Weizsacker semi-empirical mass 
formula, in which strong kinematic flavor symmetry violation is modeled by a degenerate Fermi gas. 

We check for the stability of nuclei against flssion, strong particle emission (analogous to a decay), 
and weak nucleon emission. For two light quarks with charges 2/3 and —1/3, we flnd a band of 
congeniality roughly 29 MeV wide in their mass difference, with our own world lying comfortably 
away from the edges. We also find another, less robust region of congeniality with one light, charge 
— 1/3 quark, and two heavier, approximately degenerate quarks with charges —1/3 and 2/3. No other 
assignment of light quark charges yields congenial worlds with two baryons participating in nuclei. 
We identify the region in quark-mass space where nuclei would be made from three or more baryon 
species. We discuss issues relevant to the congeniality of such worlds, but a full characterization of 
them is left for future investigation. 



I. OVERVIEW 
A. Introduction 

It is logically possible that certain physical quantities might be environmentally selected, rather than uniquely 
determined by fundamental principles. For instance, although early scientists such as Johannes Kepler sought an 
explanation for the sizes of the orbits of the planets in terms of fundamental mathematical laws, we now recognize 
that these result from complicated dynamics that involve significant environmental contingencies. Many solar systems 
are possible (indeed, many are now known to exist) with very different configurations. Today we understand that 
the size and shape of the Earth's orbit is environmentally selected: we happen to live on a planet whose orbit allows 
intelligent observers to evolve. If the orbit were much different, wc would not find ourselves inhabiting this particular 
planet. The extent to which the size of other planetary orbits in our solar system are environmentally constrained is 
a subject of active inquiry. For instance, Jupiter may play a role in stabilizing the inner solar system over large time 
scales. Other orbits, such as that of the planet Mercury, might have little to do with the emergence of life on Earth 
and consequently might not be explained by environmental selection. 

From time to time modern scientists have suggested that the values of certain physical constants, and perhaps 
even the form of some of the laws of nature, might similarly be environmentally determined P, Q. This idea, 
known loosely as the anthropic principle has received support recently on both the theoretical and phcnomenological 
fronts. Inflationary cosmology suggests that our universe may be one of a vast, perhaps infinite, number of causally 
disconnected universes, and studies of flux compactification in string theory suggest that — at least in that version of 
quantum gravity — the laws of nature and physical constants could be different in each of those universes Q . If that 
were the case, then each universe would evolve according to its own laws and only some would lead to the evolution 
of physicists who would try to make sense of what they observe. In those universes some of the laws of nature and 
physical constants may be explained not by some deep principle, but rather by the requirement that they be consistent 
with the evolution of intelligent observers 0, [1] . 
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Whatever one may think of these theoretical arguments, the prediction of the value of the cosmological constant 
A on essentially these grounds deserves to be taken seriously. At a time when the experimentally measured value 
of A was consistent with zero, Weinberg Q pointed out that values outside of a very small range (on the scale of 
fundamental physics) would prevent galactic structure formation and that probabilistic considerations suggested a 
value close to the environmental bound imposed by the structure we see around us (suggestions along similar lines had 
been made earlier by Linde 01, Banks and by Barrow and Tipler in Sec. 6.9 of 01). Shortly after Martel, Shapiro, 
and Weinberg refined this anthropic argument in a positive value of A was measured that was consistent with 
their prediction . Such an environmental explanation is all the more persuasive in the absence of any attractive 
alternative derivation of A from a dynamical principle. 

Anthropic reasoning raises some deep questions that will be difficult (perhaps even impossible) to answer scientifi- 
cally: Is there a multiverse to begin with? Are string theory and eternal inflation the correct tools with which to study 
the space of possible universes? What parameters of the Standard Model (SM) are environmentally selected, and what 
is the a priori distribution over which environmentally selected parameters range? These questions would have to be 
addressed before one could claim that some parameter, like the up-down quark mass difference, is environmentally 
selected. However, there is another, less daunting question to be asked about the parameters of the SM: To what 
extent is a particular parameter environmentally constrained? Over what range of a parameter do the laws of physics 
allow for the existence of an observer? This is the question we attempt to formulate and, at least in part, answer in 
this paper. The answer to this question might be of little scientific significance outside of an anthropic framework; 
nevertheless, the question is well formulated and can, at least in principle, be answered if we understand the relevant 
physics well enough. 

In order to determine the environmental constraints on the interactions and physical constants in our universe, we 
must be able to deduce the physical consequences of varying the laws and parameters that are familiar to us. If we do 
not understand the theoretical landscape well enough to declare certain regions environmentally barren of observers 
or others fruitful, it will not be possible to place environmental constraints on the various ingredients of the SM. The 
generally-accepted environmental explanation for the radius of the earth's orbit provides a good example of how this 
process might work: There are environmental constraints on carbon-based life forms. Biochemical structures cannot 
survive if it is too hot and biochemical reactions grind to a halt if it is too cold. One could, in principle, make these 
limits precise by studying the dependence of carbon chemistry and biology on temperature. At the end of a long 
and careful study one could imagine arriving at an understanding of the temperature constraints on the emergence 
of carbon based life forms. Combined with a similar understanding of planetary science, these would translate into 
environmental constraints on the orbital radius of a planet inhabited by carbon-based life. These questions are clearly 
within the purview of science, and this is the type of analysis we would like to carry out within the SM. The question 
of whether the earth's radius was environmentally selected depends on whether the universe contains many planetary 
systems and on the distributions of orbital radii, compositions, and so forth. Likewise, the question of whether 
some aspect of the SM is environmentally selected depends on whether there is a multiverse and on the a priori 
distribution of SM properties and parameters over that multiverse. This second question involves totally different 
physics — cosmology, grand unification, etc. — and is much harder to answer than whether some parameter of the SM 
is environmentally constrained. 

In order to carry out an environmental analysis of the SM we must first clearly specify the criterion for an an- 
thropically acceptable universe, i.e. one in which wc believe observers may exist. Then we must adopt and justify a 
specific scheme for varying ("scanning") some parameters over the landscape of possible worlds, while holding other 
parameters fixed. This procedure defines the slice through the space of SM parameters which is to be investigated. 
Then we must use our knowledge of the physics of the SM to predict the consequences of scanning those parameters 
and, in particular, to determine which choices of parameters would correspond to worlds that meet our criterion of 
anthropic acceptability.^ 

In this paper we attempt to carry out such an environmental analysis for the masses of the three lightest quarks. 
We have chosen to study quark masses because the dependence of the strong interaction on the masses of the quarks 
is non-trivial and yet is understood well enough to offer some hope that we can make definitive statements about the 
probability of the evolution of observers as a function of them. 

We propose to vary the masses of quarks holding as much as possible of the rest of SM phenomenology constant. 
In particular we leave the mass and charge of the lightest lepton, the electron, and the mass of the lightest baryon 
fiavor multiplct unchanged. In pursuit of our objective we arc forced to vary the masses of the muon and tauon (and 
other mass scales) relative to the QCD mass scale, but this has no effect on the questions of interest. We also leave 



Note that the approach wc arc taking allows only for local variations about the laws of physics that we see in our universe. We have 
nothing to say about universes with entirely different interactions and particle content. 
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the weak interactions of quarks and leptons unchanged. It is not sufficient to simply hold all the other parameters of 
the SM fixed and vary only the quark masses. 

If one insisted on varying quark masses while keeping all the other parameters of the SM fixed, for example, at 
some unification scale or at the scale of electroweak symmetry breaking, then the resulting changes to the low energy 
physics would be too extreme for us to analyze. In particular, nuclear binding in the SM depends strongly both on 
the quark masses and on Aqcd, the scale below which QCD becomes strongly coupled. If we imagined that Aqcd 
were set by the renormalization group running of a unified gauge coupling held fixed at some high mass scale Mq, 
then Aqcd would change with the masses of the quarks and leptons. 

Holding Aqcd fixed as the light quark masses are varied is not only inconvenient for calculation, but also as 
theoretically unmotivated as adjusting it in a mild way, as we prefer to do here. In Sec. |II]we describe in detail which 
parameters of the SM must vary along with the quark masses in order that we have enough control of the strong 
interactions to make sharp predictions about nuclear masses. 

The worlds we will study will look much like our own. They will have some stable baryons, some of them charged. 
The charged stable baryons, whether positively or negatively charged, will capture electrons or positrons to form 
neutral atoms^ with chemistry that will be essentially identical to ours. This leads us to adopt the existence of 
stable nuclei with charge one (some isotope of hydrogen) and charge six (some isotope of carbon) to be the criterion 
for an acceptable universe. One could argue that other elements, especially oxygen (charge eight) are necessary for 
carbon-based life. We find that worlds with stable carbon also have stable oxygen and a variety of other elements, so 
it is not necessary to refine our criteria in this regard. Note that carbon might have very different nuclear physics in 
worlds with different quark masses. If, for example, the charges of the stable baryons were different from our world, 
then the baryon number of carbon would not be 12, but the chemistry of the element with charge six would be nearly 
the same as the chemistry of carbon in our world. 

For convenience, from now on we will call a universe congenial if it seems to support the evolution of an observer. 
An uncongenial universe is the opposite. Our criterion for congeniality is, as already stated, the existence of stable 
nuclei with Z ~ 1 and Z = 6. We realize that these may not be necessary if some exotic form of life not based 
on organic molecules were possible. Likewise they may not be sufficient if there were some other obstruction to the 
development of carbon based life. The former is interesting but very difficult to study. The latter is very important 
and requires further investigation. 

For instance, Hoyle famously predicted — with better than 1% accuracy — the location of a 7.65 MeV excited state of 
^^C, based on the fact that, without such a resonance, stars would produce only a negligible amount of carbon, since 
the fusion of three a particles would proceed very slowly [llj. It is difficult, however, to translate the requirement of 
efficient carbon synthesis into a statement about the fine tuning of the fundamental SM parameters, since the rate 
of carbon production depends on the spacing between two different nuclear energy levels (which are very difficult to 
compute from first principles) and also on stellar temperatures (which could conceivably be different in other universes 
with different interactions). Furthermore, Hoyle's argument presupposes that most of the carbon available for the 
evolution of organic life has to be produced in stars by the fusion of the three a's, which might or might not be true 
of all universes in the landscape that have stable isotopes of carbon. 

Another example of an argument for an astrophysical obstruction to the development of carbon-based life comes 
from the work of Hogan, who has argued that the existence of a stable, neutral baryon would short circuit big-bang 
nucleosynthesis and/or leave behind deadly clouds of neutral baryons [l^ . However, some of the worlds we study 
have stable neutral baryons but may have different early histories than Hogan considers: for example, in some of 
our worlds the deuteron is not stable while ^H (the tritium nucleus, or "triton") is, a possibihty not contemplated 
by Hogan. The lack of a stable deuteron could interrupt the reaction path that Hogan identifies as leading to an 
explosive Big-Bang r-process in worlds with stable neutrons. 

Furthermore, the availability of stable ^H would change stellar dynamics significantly, since two tritium nuclei can 
make "^Hc plus two free neutrons by a strong interaction, in contrast to the weak interaction, pp —y^H e+ in our 
world. Finally, whether the presence, in some given universe, of roaming clouds of stable neutral baryons left over 
from Big-Bang nucleosynthesis would prevent the evolution of life everywhere presumably depends on very complicated 
astrophysical considerations. 

This highlights a distinction between two types of criteria for congeniality. The first, which we adopt, is that the 
laws of nature should be suitable for the evolution of life (as we know it) . The second, which we put aside for later 
work, is that the history of a universe evolving according to those laws should lead to intelligent life. Clearly the 
latter is a much more difficult problem, involving aspects of astrophysics, celestial dynamics, planetary physics, etc., 
each of which has its own complex rules. Returning to the planetary analogy, it is akin to showing that planets in the 



^ If the stable hadrons have negative charge, neutrinos would offset the negative lepton number of the positrons. 
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habitable zone would have any number of other attributes, like atmospheres, water, plate tectonics, etc., that might 
be necessary for life. For example, none of the astrophysics or celestial mechanics leading to earth-like planets seems 
to rely on the existence of unstable, but long-lived a-emitting nuclei. However, the earth would be a featureless ocean 
with very different life forms if it were not for plate tectonics powered by uranium and thorium decays in the earth's 
core. Such considerations are clearly important — Weinberg's constraints on the cosmological constant follow from 
consideration of the history of the universe, not the form of the SM — but they are beyond the scope of this paper. 
Constraints on the formation of life from such cosmological or astrophysical considerations can always be cleanly 
superimposed upon the constraints from nuclear physics studied in this paper. 



B. Variation of quark masses 



Quark masses in the SM, ttIq, are determined by Yukawa couplings, ga, by nia = gaV, where v is the Higgs vacuum 
expectation value (vev). In our world, the quark Yukawa couplings range from ^ 10""" for the u quark up to ~ 1 for 
the t quark. The a priori connection, if any, between the quark masses and the scale of the strong interactions, Aqcd, 
is unknown. Wc would like to consider the widest possible range of variation of quark masses, limited only by our 
ability to study their consequences with the tools of the SM. Quarks fall into two classes: light quarks, with masses 
less than or of the order of Aqcd, and heavy quarks, with masses above Aqcq. Heavy quarks do not participate in 
the physics of nuclei and atoms. If Yukawa couplings arc varied widely, the number of light quarks can range from 
zero to six, and their charges can be +2/3 or —1/3, constrained only by the fact that there are at most three of either 
charge. 

In our world, with light u and d quarks and a somewhat heavier s quark, the methods of flavor SU{3) perturbation 
theory have been moderately successful in baryon spectroscopy. Using flavor SU{S) perturbation theory wc can 
extract enough information from our world to predict the masses of the lightest baryons in worlds with up to three 
light quarks, provided their mass differences are small enough to justify the use of first order perturbation theory. 
In practice this is not an important limitation because we find that a quark species does not participate in nuclear 
physics unless its mass is considerably less than the mass of the s quark in our world. So first order perturbation 
theory in quark masses should be even better in the other worlds we explore than it is in our world. In fact, the 
principal limitation to our use of perturbation theory in the quark masses comes from possible higher order corrections 
to the extraction, from the experimentally measured baryon masses, of the reduced matrix elements of SU (3) tensor 
operators. For a similar reason, we are limited to at most three species of light quarks by the fact that we cannot 
extract flavor SU (4) reduced matrix elements from the measured baryon spectrum, given how badly that symmetry 
is broken in our world. 

To summarize, we consider worlds where baryons made of up to three species of light quarks participate in nuclear 
physics. We circumvent ambiguities associated with the extraction and renormalization scale dependence of quark 
masses by presenting results as functions of the ratio of quark masses to the sum of the three light quark masses, mr, 
in our world. 



C. Our slice through the Standard Model 

Previous research has considered the environmental constraint on the value of the Higgs vev v, with the Yukawa 
couplings fixed to the values they have in our world p^ . Along this particular slice through the SM parameter space, 
all quark and lepton masses scale together. The requirement that stable atoms exist is found to roughly impose the 
constraint 

0<^<5, (I.l) 

where v® is the Higgs vev in our universe.^ A general argument has been proposed for why this particular choice of 
a slice might correspond approximately to the landscape resulting from string theory and eternal inflation [l3 | . 

In this paper we will consider a very different slicing (see Fig. [1]), in which the three light quark masses are allowed 
to vary independently, with minimal variation on the other phenomenological aspects of the SM. First we must fix a 
scale with respect to which all other relevant masses will be defined. We do this by fixing the electron mass equal to 
its value in our world, = m® = 0.511 MeV. Next, we choose a set of quark mass values. Our world, for example, 



^ Throughout this paper we will use the super- or subscript "©" to label the value of a parameter in our own world. 
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is specified by m„ = m® = 1.5-3 MeV, rud = m® = 2.5 - 5.5 MeV, and = to® = 95 ± 25 MeV [l^. Finally, we 
adjust Aqcd so that the average mass of the lightest baryon flavor multiplct (usually an SU (2) doublet) equals the 
nucleon mass in our world, Af ® = 940 McV. 

This procedure defines our "slice" through the parameter 
space of the SM. The other parameters of the SM including 
the masses of heavy leptons and neutrinos, the parameters 
of the Cabibbo-Kobayashi-Maskawa (CKM) and Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) matrices, the scale of weak 
symmetry breaking, etc., do not impact nuclear physics sig- 
nificantly, so our analysis applies to wide families of worlds 
where those parameters vary from their values in our world. 

Because of this relatively clean separation into one set 
of mass scales that influence nuclear physics significantly 
({TOg}, TOe, Aqcd) and all the rest that do not, there is an 
additional mass rescaling that enlarges the set of worlds to 
which our analysis applies. Consider a point on our slice 
through the SM, labeled by quark masses, Aqcd, and TOg. 
Now consider a world where all these mass scales are multi- 
plied by some factor, a. Nuclear physics in this world can be 
mapped back to ours by redefining 1 MeV 1/a MeV. All 
the other mass parameters of the SM, including for example 
TO^, Tot, Mw,z, and so forth, are scaled by 1/a, but these 
have negligible effect on the question of congeniality that in- 
terests us. In effect, by fixing me we are "choosing a gauge" 
with respect to this additional invariance. 

We use AqcD to rescale the lightest baryon's mass to 940 MeV. The masses of other hadrons, in particular of the 
pion, the p(770) and the /o(600), (or the equivalent mesons in worlds with different light quarks), which determine the 
properties of nuclei, are affected by the choice of quark masses and of Aqcd- Changing these mesons' masses affects 
the binding of nuclei and introduces uncertainty into our results. As a result there are some regions of quark mass 
parameter space where we can make reliable quantitative statements and others where only qualitative statements 
are possible. This question is discussed in considerably more detail in Sec. Ill Dl 

We have chosen a slice through the parameter space of the SM along which it is relatively simple to examine 
congeniality as the quark masses vary. This is only one of many possible slices, and might not be the slice dictated by 
a deeper understanding of the landscape of possible universes. We choose our approach simply because the resulting 
universes can be studied with some certainty. So doing, we accomplish two things: first, we give an example of what 
is required to analyze the environmental impact of non-trivial variations in SM parameters, and second, we find some 
of the limits of the domain of congeniality in which we live, and also discover some different (disconnected) domains 
of congeniality elsewhere in light quark mass parameter space. As our knowledge of QCD improves, in particular 
as lattice QCD matures, it may become possible to explore much larger domains of parameter space by simulating 
hadrons and nuclei on the lattice. 



D. Plan of the paper and survey of results 

In Sec. HI] we confront some of the (meta)physical issues that arise when attempting an environmental analysis 
of quark masses. First we define more carefully the difference between environmental constraint and environmental 
selection. We discuss what features the a priori distribution of quark masses might have. Next we explain why we 
ignore quarks with masses greater than the confinement scale of QCD. This leads us to a parameterization of the 
space of light quark masses and finally to some discussion of baryon masses, nuclear forces, and stability. 

In general, only members of the familiar fiavor SU (3) octet of spin-1/2 baryons can play a role in nuclear physics.^ In 
Sec. Illll we express the masses of the octet baryons as functions of the three light quark masses. Most of the information 
necessary to relate octet baryon masses to quark masses can be obtained from a study of octet baryon mass differences 
in our world, using first order perturbation theory in the quark masses. However, baryon mass differences cannot tell 




FIG. 1: Two slices through the parameter space of the SM are 
shown in this cartoon representation of that space. For sim- 
plicity the space is drawn as two-dimensional and the slices as 
one-dimonsional. The two slices meet at the point the corre- 
sponds to the parameters of our own world (labelled "us" in 
the figure). We may imagine slice 1 to correspond to varying 
quark masses at fixed as{Mz), while slice 2 would correspond 
to the subspace we investigate in this paper. 



* Two exceptions are mentioned: when all quark masses are small one must consider the possibility that the dihyperon is the most stable 
baryon (see Sec. IV ETI i. and when only one quark is light, the spin-3/2 decuplet must be considered (see Sec. IV At . 
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us about the dependence of baryon masses on the average quark mass. We take this information from the tr-term in 
ttN scattering (see [l^l and references therein) . To get as accurate an estimate as possible of the connection between 
baryon masses and quark masses, we include estimates of the electromagnetic contributions to baryon masses. 

Except in the case where all three light quark masses are nearly equal, only the two lightest baryons participate in 
nuclear dynamics. The others, like hyperons in our world, make only exotic, very short lived nuclei (usually called 
"hypcrnuclei" ) When all three quark masses are nearly equal, then all eight members of the baryon octet may 
be involved in nuclear dynamics. In Sec. IIVI we estimate the masses of nuclei in our slice through the landscape 
as functions of the baryon masses. For some universes a given nuclear mass can be estimated very accurately by 
considering its analog: a nucleus in our world whose binding energy should be the same except for easily calculable 
Coulomb corrections that occur if the participating baryon charges differ from those in our own world p^ . For 
other universes where direct analogs are not available, we build a generalized Weizsacker semi-empirical mass formula 
(SEMF), with parameters more or less closely related to those in our world. Kinematic flavor-breaking contributions 
are obtained from a degenerate Fermi gas model [20| . supplemented by a term proportional to the quadratic Casimir 
of flavor SU{2) or S'C/(3), depending on whether there are two or three very Hght quarks. 

With these tools in hand, in Sec. |V]we explore the congeniality of worlds along our slice through the SM. Since 
nuclear masses can be estimated, it is possible to check for the stability of various nuclei against fission, a decay, and 
"weak nucleon emission" [i.e., the emission of a nuclcon and a lepton-antilcpton pair, mediated by a weak interaction). 
For universes with only two light quarks and charges -1-2/3 and —1/3, u and d, we find a band of congeniality roughly 
29 MeV wide in m„ — TOcj, with our own world lying comfortably away from the edges. We also find another region of 
congeniality with small rud and m„ « TOs (or equivalently, small and m„ sa rud), whose exact width in to„ — rus 
is difficult to estimate because nuclear forces in these worlds differ qualitatively from ours. Worlds with two light 
quarks of the same charge^ either -|-2/3 or —1/3, are not congenial since, regardless of the quark masses, carbon is 
always unstable. 

Finally, we identify the region in the space of quark masses where nuclei would be made from three or more 
baryon species. We discuss some issues relevant to nuclear stability in such worlds, but a full characterization of 
them as congenial or uncongenial is left for future study. We point out several possible sources of uncongeniality for 
these worlds: it is possible that in these worlds octet baryons would decay into dihyperons, which are bosons and 
would therefore behave very differently from micleons in our universe. In worlds with one positively-charged and two 
negatively-charged light quarks, nuclei with net electric charge would be very heavy, and might well be unstable to 
decay by emission of light, uncharged nuclei composed of a pair of each baryon species. 

II. SETTING UP THE PROBLEM 

In this section we delve more deeply into some of the issues raised in Sec. [H We argue that a study of the 
environmental constraints on quark masses is a well-posed problem, and explain its relation to more ambitious 
anthropic inquiry. Also the slice through the SM parameter space introduced qualitatively in Section [T] must be 
delineated more quantitatively. 

A. A priori distribution of quark masses 

As explained in Sec. [H we focus on the relatively well-defined problem of environmental constraints on quark 
masses. Nevertheless we recognize that there is much interest in the theoretical physics community in whether physical 
parameters might be environmentally selected. The jump from constraint to selection requires some theoretical input 
(or prejudice) about the a priori distribution of quark masses. Also, a prejudice about the a priori quark mass 
distribution can suggest where to focus our analysis of environmental constraints. Therefore, we review here the 
relation between environmental constraint and environmental selection, and then summarize our thoughts about 
possible a priori quark mass distributions. 

Anthropic reasoning can be presented in the language of contingent probabilities [2l| . According to Baycs's theorem, 
the probability distribution for measuring some set of values for the quark masses can be expressed as 

p({mi} I observer) cx p(obscrvcr|{mi}) x p{{mi\) , (H-l) 

where p{{mi}) is the a priori distribution of quark masses over the landscape. The quantity p({TOi} [observer) is the 
contingent probability distribution for the quark masses over the subset of the landscape that contains intelligent 
beings likes us, capable of measuring quark masses. In simple terms, it is the probability that some set of quark 
masses would be measured. The quantity p(observer|{TOi}), the contingent probability for an observer existing, given 
a choice of quark masses, is the stuff of environmental constraint and the primary subject of this paper. 
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FIG. 2: The quark masses and AqcDi shown on a logarithmie scale. 



The connection between environmental constraint and environmental selection is provided by the a priori distri- 
bution, p{{mi}). Environmental selection requires the existence of an ensemble of universes over which p{{mi}) is 
defined. Anthopic reasoning is based on some prejudice about the form of p{{mi\) as (presumably) determined by 
some fundamental theory, such as string theory in an eternally-inflating scenario. This should be contrasted with a sit- 
uation where one or more quark mass is dynamically determined. Then p{ma) (for some a) would be a delta-function 
and ma would not be environmentally selected. 

In his analysis of the cosmological constant A, Weinberg argued that the a priori probability distribution p(A) 
should not vary significantly over the range where the contingent probability p(obscrvcr|A) is non-zero, in which case 
the a priori distribution factors out of the expression for p(A| observer). Similarly, if p({mi}) does not vary significantly 
over the region where p(observer|{?7ii}) is non-zero, then the probability that we observe some set of quark masses, 
which is the object of fundamental interest in an anthropic context, is simply proportional to p(observer|{mi}), which 
we analyze here. 

We believe that if anthropic considerations have any relevance at all, then there is reason to assume that the 
logarithms of the quark masses are smoothly distributed over a range of masses small compared to the Planck scale. 
Note, first, that an environmental analysis is predicated on the absence of any other underlying explanation for quark 
masses. Therefore we assume that the probability of a pattern of masses is the product of independent probabilities, 



In this paper we will usually [i.e., unless otherwise specified) label the positively-charged quarks as u, c, and t, and 
the negatively-charged quarks as d, s, and &, in ascending order of mass. (Therefore the identification of the label 
1 = 1, 6 in eq. pi.2[) with u, d, s, c, b, t varies according to the values of the {mi}.) In general, we shall care little 
about the details of the fiavor structure of the weak interactions, but in Sees. IIVI and |V] we will assume that no 
element of the CKM matrix is accidentally zero, so that no baryon that would normally decay by weak interactions 
is accidentally stable. 

Taking an anthropic point of view, there is a hint that, at least at masses below the scale of electroweak symmetry 
violation, the a priori distribution of the quark masses might be logarithmic. This hint comes from the observed 
distribution of quark masses in our world, which is shown in Fig. [21 The s, c, b, and t quarks are so heavy (on the 
scale set by Aqcd) and short lived that they seem to play no role in nuclear and atomic physics (and consequently 
not in chemistry or biology either). Therefore, referring to Eq. pi.ip . it is reasonable to take p(observer|{mi}) to be 
essentially independent of m^, TOc, fnt, and mt. Then the pattern of heavy quark masses that we see could well be a 
measure of the a priori probability distribution p{m). 

Note, though, that the authors of [l^] proposed that m* might be anthropically selected to be near its measured 
value, ~ 170 GeV. This proposal was based on the assumption (unmotivated by any dynamical consideration about 
the landscape) that the a priori distribution for the Higgs quartic self-coupling, A, strongly favors small values. The 
environmental requirement that the SM have a stable or mctastable vacuum would then select mt through its effect 
on the renormalization group running of A. If we relax the assumption about the a priori value for A, however, there 
would not seem to be a reason to expect strong anthropic constraints on mt. In any event, such a constraint could 
only apply to the heaviest quark, leaving the masses of s, c, and b unconstrained. 

Although we have only four (or three) data points, the measured values of heavy quark masses range over more 
than two orders of magnitude, which suggests that the a priori cumulative distribution 
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(11.2) 




(II.3) 



might be a smooth function in the logarithm of the mass, over some range that contains the six quark masses measured 
in our world [23| . The a priori distribution of quark masses is determined by the landscape, which presumably is 
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generated by physics at the Planck scale. Therefore the scale for variation in the quark mass probability distribution 
should be Mpjanck and we expect 

771 

P(m) = a log- + ... (II.4) 

-iWpianck 

where the omitted terms vary more slowly with m than the logarithm. 

Keeping only the leading term in Eq. pi.4p implies that quark masses in that range are a priori equally likely to 
fall between any two consecutive orders of magnitude, and corresponds to a scale-invariant distribution 

p(rn) oc — . 

771 

Approximately scale invariant distributions are observed in a variety of phenomena associated with complex underlying 
dynamics: e.g., the distribution by area of lakes in Minnesota has such a form, with a lower cutoff given by the 
conventional definition of a lake and the upper cutoff given by the area of Minnesota.^ In the context of a possible 
landscape of quark and lepton masses, a scale-independent a priori distribution was investigated in [2^. 

Note that an a priori logarithmic distribution of quark masses puts a much greater weight on masses that are very 
small compared to Aqcd than would a linear distribution. Therefore it is interesting to explore the congeniality of 
worlds where all light quark masses are very small compared to Aqcd- We have been able to make some progress in 
formulating the tools needed to explore this region of parameter space, and we can make some qualitative statements 
about its congeniality. Nevertheless we have not been able to explore this region as fully as we would like. 

Finally, we emphasize again that this excursion into a discussion of a priori probabilities is independent of the 
analysis of environmental constraints on quark masses which will concern us in the rest of this paper. For a more 
explicit model of a possible landscape of Yukawa couplings, coming from the overlap of localized wavefunctions 
randomly distributed over compact extra dimensions, see [26|. 



B. Light and heavy quarks 

The energy scale at which QCD becomes strongly coupled, Aqcd, provides us with a distinction between light and 
heavy quarks. Heavy quarks are those with rrij > Aqcd and light quarks have mj < Aqcd- Even though the strange 
quark mass in our world is heavy enough that it decouples from nuclear physics and chemistry, it is nevertheless light 
compared with Aqqj-,, so it falls among the "light" quarks in this discussion. If any quarks are lighter than Aqcd, 
then assuming no accidental zeros in the CKM matrix, the heavy quarks will decay rapidly (on the time scales of 
atomic physics or chemistry) into light quarks by means of weak interactions. In this case only light quark masses 
would be environmentally constrained by nuclear physics, atomic physics, and chemistry — the constraints explored in 
this paper. 

If all quarks were heavy, then the emergence of an observer would require a dynamical explanation, which would 
be an interesting counterexample to anthropic reasoning: an observer who found him/her self in a world where all 
quark masses are large compared to Aqcd would have to dismiss an anthropic origin and instead seek a dynamical 
explanation for the pattern of quark masses. 

To see why, suppose c were the lightest quark in a world with all quark masses greater than Aqcd - If no other quark 
had mass close to nic then only the ccc baryon (analogous to the A"*""*" for Qc = 2/3 or to the fl^ for Qc = —1/3) 
would participate in building nuclei. It is easy to see that Coulomb and symmetry effects would destabilize all systems 
with baryon number greater than one (remember, no single species nucleus, starting with the diproton or dineutron, 
is stable in our world), and there would be no congenial worlds. 

Suppose, on the other hand, that two heavy quarks had masses close enough to one another that both could 
participate in nucleus building. The most favorable situation, from the point of view of Coulomb interactions, occurs 
when the quarks have different charges, so to make the discussion concrete we can call these heavy, nearly degenerate 
quarks the c and b quarks. The two lightest baryons, ccb and bbc, analogs of the proton and neutron, would experience 
attractive forces with a range ~ 1/Aqcd-^- Attractive forces between massive particles bind when the range of the 



^ Scale invariance also explains the well-known phenomenon — known as Benford's law — that the first digits of numbers from many different 
kinds of real-life data are logarithmically distributed |24|| . 

® Although there are no mesons with masses of order Aqcjj in such a world — the lightest meson would be fee, the analog of the — 
nevertheless the nucleons would have radii of order Aqcd (up to logarithms), and when they overlap, they would attract by quark 
exchange processes. This apparent violation of the rules of effective field theory (the longest range force is not associated with the 
nearest threshold in the complex plane) was previously encountered in heavy meson effective field theory and is explained in Ref. |27|I . 
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(a) (b) 



FIG. 3: (a). Graphical representation of the landscape of light quark masses for fixed rriT = rau + 'tIjj + ms; (b). Reduced landscape, 
assuming rris > m^. 

force is much greater than the particle's mass''. So it is reasonable to expect stable nuclei when m;, « rric > Aqcd- 
However the two heavy quarks' masses would have to be within a few tens of MeV of one another, otherwise the nuclei 
would rapidly decay away by weak nucleon emission. This process is discussed in Sec. IIVB 51 and the corresponding 
constraints on the baryon masses are covered in Sec. IV CI For now, it is sufficient to know that the quark mass 
difference cannot be greater than of order tens of MeV. 

Thus there is a possible window of congeniality if two heavy quarks are nearly degenerate. However, if the a priori 
quark mass distributions, p{mj) are independent, then this is a subspace of very small measure in the ensemble of 
universes. If p(m) is a smooth function of logm up to the Planck scale, then this congenial domain has essentially 
measure zero. Physicists who lived in that universe would be confronted by an "accidental" symmetry, similar to (but 
even more exaggerated than) the one that puzzled particle physicists in our universe before the discovery of QCD. 
Before circa 1970, it was thought that isospin was a good symmetry because the u and d quarks, both believed to 
have masses of order 300 MeV, had a mass difference of only a few MeV. This approximate degeneracy appeared 
mysterious, and physicists sought a dynamical explanation. When QCD explained that the mass scale of hadronic 
bound states was set by Aqcd, not by u and d quark masses, then the mysterious approximate equality was replaced 
by the less mysterious inequality, mu,md <C Aqcd- 

Analogously to the pre-1970 situation in our universe, physicists in the heavy quark world would seek a dynamical 
reason for the approximate equality of heavy quark masses. In environmental terms, they would have to explain why 
the a priori quark mass distribution should have a (approximate) S function, S(mc — rub), since a world with only 
heavy quarks would necessarily be uncongenial, unless there were a nearly-perfect degeneracy between the masses of 
the quarks, and such a degeneracy would not admit a compelling environmental explanation. * 

C. Parameterizing the masses of light quarks 

As mentioned in Subsection II CI we will consider the variation of three light quark masses over a slice through the 
space of parameters of the SM in which the average mass of the two lightest octet baryons is fixed with respect to 
the mass of the electron. 

To be definite we assume the charge assignments (2/3,-1/3,-1/3). The extension to the other mixed charge 
possibility, (2/3,2/3,-1/3), is mentioned as needed. We denote the quarks {mu,rnd,ms) in the former case and 
{mu,mc, rud) in the latter case. The case where all three light quarks have the same charge, either +2/3 or —1/3 docs 
not yield a congenial world with one possible exception: If all three quarks have charge —1/3 and all quark masses 
are much lighter than Aqcd, then it is possible that such a world would be congenial. We mention this regime in 
Section El 



^ The criterion is mVoR^ > constant, where Vb and R are characteristic of QCD and independent of m (up to logarithms), while m is 
assumed to be large. 

* Of course, just as we do not know how small quark masses can be in the landscape, it is possible that physicists in a universe with only 
heavy quarks would not realize that quark masses are dynamically allowed to fall below Aqcd- Until they had learned, as we already 
know, that the a priori distribution allows quarks to be lighter than Aqcd, they might accept the anthropic argument that only nearly 
degenerate (heavy) quark masses lead to a congenial world. 



10 




d niT 



For a fixed value of tut = ruu + + ms, we may rep- 
resent the landscape of light quark masses by the points in 
the interior of an equilateral triangle with altitude of vtlt, as 
shown in Fig. EJa). The value of each quark mass is given 
by the perpendicular distance from the point to the corre- 
sponding side of the triangle. (Recall that the sum of those 
perpendicular distances from any point inside an equilateral 
triangle is always equal to the altitude of the triangle.) 

As we first discussed in Sec. Ill Al we shall not need to make 
any assumptions about the structure of the CKM matrix 
beyond it not having accidental zeros. This removes any 
distinction between the worlds described by points in the 
right and left halves (related by d s) of the triangle of 
Fig.jS^a), so we will define the s quark to be the more massive 
of the two charge -1 /3 quarks and (when appropriate) restrict 
our attention to the reduced landscape shown in Fig. [3Kb). 
The triangular landscapes relevant for quark charges (2/3,2/3 
"T-c, fTid TTi-u, TTiu with thc convcution that nic > 

The full space of light quarks is, of course, three-dimensional, and may be represented as a triangular prism, as 
in Fig. m The altitude of any triangular slice is rriT, the sum of the three light quark masses. Remember that the 
Aqcd has been adjusted so that the average mass of the lightest flavor multiplet of baryons is held fixed at 940 MeV. 
In practice this multiplet is always a flavor doublet. The prism shown in Fig. [4] summarizes the range of universes 
explored in this paper. 

The variation of baryon masses and nuclear forces over the prism is influenced by the adjustments of Aqcd- Those 
adjustments arc described in the following subsection, and the resulting variation of baryon masses is discussed in 
Sec. 



FIG. 4: Representation of thc full three-dimensional space of 
light quark masses as a triangular prism, with four different 
slices of constant rriT = ruu + md + ms = -v/G'tio shown. 



-1/3) can be obtained by the replacements: 



D. Strength of the nuclear interaction 

As described in Sec. HI our aim is to explore universes with at least one and at most three light quarks, in which 
nuclear physics resembles ours as closely as possible. In this section we discuss more carefully the extent to which 
this can be accomplished. 

We find three regimes where long lived nuclei can form: 

(i) Universes where all three quarks are very light compared to Aqj^^; 

(ii) universes where one quark is much lighter than the other two, which are approximately degenerate; and 

(iii) universes similar to ours where two quarks arc considerably lighter than thc third. 

The regions are sketched on the quark mass triangle in Fig. [5l These regions are only well separated when quark 
masses are large. When all quark masses are small compared to the ~ 10 MeV scale of nuclear binding, then only 
case (i) can occur. This is shown in Fig. [Ub). Other universes outside of regions (i), (ii), or (iii), for example those 
with a hierarchy like mi <C m2{^ 10 — 20 McV) <C ^3 < Aq^j-, have no stable nuclei. 

In universes of type (i) all eight pseudoscalar Goldstone bosons (tt, K, and r]) participate in nuclear interactions. 
This makes it impossible to describe nuclei quantitatively in those worlds based on information from ours. In universes 
of type (ii) there is an approximate flavor SU{2) symmetry, but unlike our world, the lightest pseudoscalar Goldstone 
boson is a flavor singlet and there is a complex doublet analogous to the K/K system that is not much heavier. Again, 
nuclear physics in these worlds will differ quantitatively from ours. It is still possible to make qualitative statements 
about nuclear physics in these regions because the qualitative features of baryon-baryon forces should not change: 
they are short range, attractive at long distances and strongly repulsive at short distances. 

In the last case, (iii), the lightest pseudoscalar bosons form a triplet under the SU{2) symmetry associated with 
the two lightest quarks. In this very important case, nuclear physics is quantitatively similar to our world over a 
significant range of quark masses. Without loss of generality, we can call the two lightest quarks u and d and the 
heavier one s. The rest of this subsection is given over to a discussion of case (iii). 

Nuclear forces depend on the masses and coupling constants of the mesons that couple strongly to the light baryons. 
In case (iii) this includes the pion, which gives the long range tail of the attractive nucleon-nucleon force, the /o(600), 
also known as the a, which is a short hand for the correlated two pion exchange forces that generate the intermediate 
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FIG. 5: (a). Three regions in quark mass parameter space in which complex nuclei form. The sum of quark masses here is large compared 
to the energy scale of nuclear binding. The red shaded regions are not congenial because both light baryon species have the same charge 
(see Sec. [Vj. (b). In a slice where the sum of quark masses is smaller, the peripheral regions disappear, leaving only the central region 
where all quark masses are small compared to nuclear binding energies. 



range attraction in the nucleon-nucleon force [28[, and the p(770) and w(783), which generate the short range repulsion. 
The same changes of parameters that keep the mass of the nucleon fixed over regions (iii) of the quark mass prism, 
tend also to keep the masses of hadrons that scale linearly with quark masses fixed as well, thereby reducing the effect 
of quark mass variations on the nature of the nuclear force. 

If it were not for the pion, we would argue that non-strange meson masses and coupling constants are so insensitive 
to the small variations in the u and d quark masses, that the variation of the nuclear force could be ignored over the 
entirety of type (iii) regions. The pion mass, however, varies like the square root of the average u and d quark masses, 
rh ~ {rriu + md)/2, and therefore changes more significantly over the prism. For example, as to ^ the pion masses 
get very small, although the p, ui, and /o(600) masses hardly change at all (recall that we are adjusting Aqcd to keep 
Mat fixed as we send to ^ 0.) 

Let us make this more precise: We are only interested in quark masses that are small enough compared to Aqcd 
so that hadron masses, other than those of the pscudoscalar Goldstone bosons, can be accurately estimated by first 
order perturbation theory. 

Ma ~ OqAqcd + bam + Carris (II. 6) 

for hadron a. Here Oq,, ba, and Ca are dynamically determined constants that arc independent of quark masses. The 
triplet of pscudoscalar bosons, which we will call pions, are different. Their masses vary dramatically for small quark 
masses, as dictated by the Gell-Mann-Oakes-Renner relation [29j.^ 



= , (11.7) 

where the pion decay constant /^r and the chiral-symmetry breaking vev of QCD (qq) are believed to be nearly 
independent of light quark masses. So schematically, 



rrirr = rfTTV^AqcD rh ■ (11.8) 
It is convenient to express the hadron and quark masses in terms of their values in our world, 

VuA,s = muA,s/m®^^,., y = m/m®, yT = mT/m^ 

yN = MN/M®, = Aqcd/A®cd> y^ = m^/m®. (II.9) 



Then 



TO^ = VyA^ m® , (11.10) 
Mn = (783 MeV) y^ + (35 McV) y + (124 MeV) , (11.11) 
™p2/A + + ^IVs ■ (11.12) 



TOn 



The mass of the /9(770), a representative light meson, is decomposed in a similar way, with unknown coefficients. We 
have borrowed the numerical coefficients in Eq. (jll.lip from the fit performed in Sec. IIIII 



In this discussion we ignore electromagnetic mass differences. 
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The numerical coefficients in Eq. pLll|) reflect the smaU contribution to the nuclcon mass from u and d quarks in 
our world (35 MeV), the rather surprisingly large — though very uncertain — contribution from the s quark (124 MeV), 
and the dominant QCD contribution (783 MeV). We are interested in how meson masses vary for values of rh and rus 
that remain in region (iii), where m rris. As rh moves away from its value in our world, Aqcd must be adjusted to 
keep the nucleon's mass constant at 940 MeV. The masses of non-chiral mesons like the p{770) and /o(600) can be 
expected to vary rather little as m changes because their proportional dependence on Aqcd and rh should be similar 
to the nuclcons'. 

The pion mass scales differently with m and Aqcd, and it cannot be expected to remain nearly constant as m varies. 
So the extent to which nuclear physics changes over the quark mass parameter space of interest (type (iii) regions in 
the prism) can be examined in two stages: (a) how much does the pion mass change? and (b) how important is the 
mass of the pion in determining the parameters of nuclear binding? The first question can be answered quantitatively, 
but the second remains a matter of speculation. 

To answer the first question it is helpful to combine the Eqs. pi.lO[) and (jll.ll[) . 

Mn = (783 MeV) ^ + (35 MeV) y + (124 MeV) ■ (11.13) 

y 

It is clear from this equation that there is a subspace of the quark mass parameter space on which both the nuclcon and 
the pion mass remain fixed as rh and rriT vary. To display this constraint more usefully, we set M^v = = 940 MeV, 
and substitute values of quark masses from our world, where rh® « 3.6 MeV and mf « 93 MeV^^ into Eq. (|IlT3|) . 
and switch variables to rh and rriT, 

940 MeV = + 7.0„"., + 1.3... , (11.14) 

m 

where rh and tot are in MeV. We leave the dependence in Eq. pi.l4[) explicit so we can study the constraint 
between rh and rriT as TOtt varies. 

Equation pi.l4p constrains the quark masses to a two dimensional subspace of the prism. Since the pion mass 
depends only on rh, the constraint is independent of the difference of the light quark masses, and defines a line parallel 
to one side in a fixed nir triangle, as shown in Fig. ^a). The subspace is also constrained by the demand that the 
third (in this case, the s) quark must be heavy enough so that strange baryons decouple from nuclear physics. In 
Fig. [6jb) we plot the average light quark mass, to, that satisfies Eq. pI.14p as a function of tot- The value of rh 
does not depart far from its value in our world because a small variation in to compensates for the variation in tot in 
Eq. pi.l4[) . The shaded band in the figure shows the range of rh over which TO,r varies by a factor of two. The plot 
is truncated on the left by the constraint that the mass of the strange quark must be greater than to by ~ 5 MeV, 
to insure that the domain remains a type (iii) region. (The choice of a 5 McV splitting is justified by Eq. (jV.20p . in 
Sec. Eel) 

If nuclear binding were dominated by singlc-pion exchange, Eq. pi.7p would imply quite a strong dependence of the 
nuclear interaction on the masses of the light quarks. In particular, increasing the masses of the u and d quarks in 
our universe would be expected to significantly decrease and eventually turn off nuclear binding 0, [l^. It is believed, 
however, that the attractive forces responsible for nuclear binding receive significant contributions from correlated 
two-pion exchange in the region of 400 - 600 MeV invariant mass, an effect that is usually associated with the /o(600) 
meson [28j . The mass of the /o(600) is not singular in the chiral limit, nor docs it vary dramatically when the pion 
mass is increased well above its value in our world [lOl , so there is reason to believe that nuclear binding depends less 
dramatically on light quark masses than the single pion exchange mechanism would lead one to believe. 

Most studies of the effect on the nuclear interaction of changing the pion mass have focused on the S-wave nucleon- 
nucleon scattering lengths and on the related problem of dcuteron binding. Various attempts have been made to 
compute these numerically, either by solving the Schrodingcr equation with approximate potentials (3l] . [33 | or from 
lattice QCD calculations [33|. The exact dependence of deuteron binding on the quark masses, however, still remains 
uncertain. 



^" The errors on the quark masses are very large, but since this analysis is only intended to yield qualitative conclusions we will ignore 
them in this section. 

There are two solutions of Eq. I|II.14|| for rri as a function of j/^ and rriT- Wc can identify the relevant one by requiring m® = 3.6 McV 
for m® = 100 MeV. 

Another approach to studying the dependence of nuclear binding on the pion mass is to consider the analytic structure of nucloon- 
nucleon scattering amplitudes and to study the effects — on the corresponding dispersion integrals — of varying the pion mass l34l Issll . 
This approach might be subtle because, as was explored in l27t in a different context, the expectation from effective field theory — that 
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FIG. 6: (a) The quark mass triangle for rriT = 40 MeV showing the the line along which m-K = m® = 140 MeV with shading 
corresponding to pion mass variation by a factor of two above and below its value in our world. Note that the entire region 
satisfies rus ^ rh. (b) The constraint on m and rus as functions of rriT, generated by fixing and Mjv to their values in our 
world. As before, the shading correspond to varying by a factor of 2. The shaded regions are truncated when rUs — rh < 5 
MeV, because the corresponding worlds would have nuclei containing strange baryons. 



In our world, the deuteron plays an important role in primordial nucleosynthesis, but it is unclear whether a 
deuteron bound state is necessary in all worlds to ensure that heavy elements eventually become abundant enough 
to support intelligent life (as we have noted in Sec. II Al tritium might serve as an acceptable form hydrogen and its 
binding is more robust than the deuteron's). For this reason, an in keeping with the philosophy of this investigation 
as outlined in Sec. II A| we will not focus specifically on the deuteron binding, but rather with the stability of any 
isotopes of hydrogen, carbon, and oxygen. 

Based on general considerations of dimensional analysis, we expect that all of the mass scales Mi relevant to nuclear 
physics — including the masses of the baryons as well as the masses of the /o(600), p(770), and w(783) resonances — 
should vary roughly as 

A®cD 

where Am^ is the quark mass variation being considered. ^'^ The impact of the baryon mass variations on nuclear 
structure, however, will be much enhanced compared to the eiTcct of the variation of other mass scales. In fact, if 
it were not for the sensitivity of nuclear structure to the baryon masses, the breaking of SU{2) and SU{3) flavor 
symmetries by the strong interaction would be largely irrelevant to a nuclear physicist. 

As an illustration of the exquisite sensitivity of nuclear structure to the baryon masses, consider the fact that (as 
we shall work out in detail in Sec. IV C 4[) the A baryon could be as little as 20 MeV heavier than the proton and the 
neutron and still not participate in stable nuclei. Thus, a breaking of flavor SU (3) by a mere 2% in the baryon masses 
would leave almost no trace of that symmetry in the composition of stable nuclei. It is therefore lucky that baryon 
masses can be reliably estimated as functions of the light quark masses using first-order perturbation theory in fiavor 
SU{3) breaking (as we shall work out in Sec. IIIip and that the sensitivity of nuclear structure to baryon masses can 
be qualitatively understood by modeling the nucleus as a free Fermi gas subject to a confining pressure (as we shall 
explain in Sec. IIVB 2p . Thus, in this paper we shall retain the effect of quark mass variations on the baryon masses, 
but we shall mostly assume that we can ignore the variation in the strength of the nuclear interactions due to the 
relatively small changes to the quark masses that we consider. 



low-energy scattering is controlled by the nearest singularity in the amplitude (and therefore, in this case, by the pion mass) — can fail. 
For the 5- wave channels, the effective nucleon-nucleon contact interactions at low energies are not of order 1/m^, as naive effective 
field theory might suggest, but instead are significantly larger (as can be seen from the size of the absolute values of the corresponding 
scattering lengths). The presence (or absence) of a spin-1 nucloon-nuclcon bound state (the deuteron) could, through Reggc pole theory, 
alter the asymptotic behavior of the dispersion integrals that connect the low-cncrgy scattering amplitude to the poles and branch points 
set by the pion mass [3^ . 

In fact, investigations by [37| and [38l | suggest that, for the mass of the p(770), non-linearities in Am^ ^ are relevant near the mu = = 
chiral symmetry point, but these non-linearities actually make the quark mass dependence of Mp milder than Eq. I III.151 would predict. 
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III. OCTET BARYON MASSES 



111 this section we use flavor and chiral symmetry violation among the baryons in our world to find the masses of 
the light baryons as functions of the quark masses. These are used as inputs in Sec. llVl to determine the spectrum of 
nuclear masses, which in turn allow us to study the congeniality of worlds with different quark masses. 

In our world the lightest baryons are in the octet representation of the flavor SU (3) symmetry group associated 
with the u, d, and s quarks. These are the nucleons (proton and neutron), the E's, the S's, and the A. Quark mass 
effects (which break flavor SU{3)) and the electromagnetic interaction give rise to mass differences within the octet. 
Electromagnetic contributions to the masses are small and calculable [3^. They can be subtracted from the total 
baryon mass to leave only the contributions from QCD and from quark mass differences. 

To first order, the quark mass contributions are proportional to the reduced matrix elements of SU{i) tensor 
operators. These matrix elements can be extracted from the measured mass differences in our world. These can then 
be used to compute the masses of the octet baryons in worlds with quark masses that differ from those in our world. 
One complication is that the analysis of baryon mass differences does not enable us to estimate the matrix element 
of the flavor singlet operator. This, however, can be extracted from the analysis of chiral symmetry violation in the 
baryon sector, through the well-known cr-term analysis. 

As we pointed out in Sec. II Bl SU{3) flavor symmetry will be only very mildly broken in most of the worlds 
studied in this paper, and first-order perturbation theory will therefore be an excellent approximation. Higher-order 
corrections principally affect the extraction of the SU{S) reduced matrix elements from the measured baryon masses. 
However, we shall neglect such corrections in light of the level of accuracy that is called for in this kind of investigation. 



A. Flavor SU{3) and baryon masses 



The mass of a baryon is given by the eigenvalue of a Hamil- 
tonian H, which we may decompose as 

H = Hqcu + -ffflavor + ^^EM , (IH-l) 

where Hqcb corresponds to QCD interactions and is inde- 
pendent of quark masses and electric charges, i?flavor is the 
quark mass dependent piece 
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and Hem gives the electromagnetic contribution to mass of 
the baryon. Hem is the only part of H that depends on the 
electric charges of the component quarks. Using the renor- 
malized Cottingham formula [i^, |4l|, it is possible to re- 
late baryon matrix elements of Hem to experimental data 
for electron-baryon scattering. Data is only available for 
electron-nucleon scattering, but there is evidence that (Hem) 
is dominated by the Born terms which in turn are dominated 
by the charge and magnetic moment contributions, which 
are known. Here we use the values as quoted in Gasser and 
Leutwyler in They are tabulated in Table HI 

The quark mass dependent term in H can be decomposed in terms of irreducible tensor operators 



FIG. 7: The SU{3) weight diagram for the octet baryons 
showing their conventional names and their assignment of hy- 
perchargc Y and the third component of isospin I3. The two 
central states, denoted Chigh and Ciowi ™ix at lowest order in 
5(7(3) perturbation theory. 



-ff flavor = ITlo'dl 



(in.3) 



with 



6o = Y - {uu + dd + ssj , 

Too = —j={mu +md + nis), 
V6 



O3 = UM — dd, 



e| = Y 2 {uu + dd-2s 



2V3 



(to„ + nid - 2ms) 



(IH.4) 



Ref. [I7II does not list an electromagnetic contribution to the A mass. We have computed it using their methods. 
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Baryon 
species 



Flavor-dependent 
contribution 



Electromagnetic 
correction (MeV) 



Experimental 
(MeV) 



. ittet 
MeV 



Residual 
(MeV) 



E+ 
E" 

Chigh 
Clow 




,ms + (F + D) 7713 
=^^-°^m8-(F + D)m3 
mg + (F - D) ms 
I ms - {F - D) ms 
mg + 2F 7713 
mg — 2F 7713 



+0.63 
-0.13 
-0.07 
+0.79 
+0.70 
+0.87 
-0.21 
+0.21 



938.27 
939.57 
1314.83 
1321.31 
1189.37 
1197.45 
1192.64 
1115.68 



939.87 
942.02 
1316.81 
1323.38 
1188.42 
1197.21 
1191.82 
1109.61 



1.60 
2.45 
1.98 
2.07 

-0.96 
-0.24 
-0.82 
-6.07 



TABLE I: Contributions to the octet baryon masses in our world. The flavor-independent contribution is Ao = (Hqcu) + 
7710 (©o)- The quark mass parameters mo,3,g are defined in Eq. HIII.4[I . 



The Hamiltonian, Eq. (jIII.3p violates SU{3) symmetry but commutes with the generators Q3 and Qs-^^ Therefore, 
no matter how light quark masses vary, we may label octet baryons by the eigenvalues of the third component 
of isospin, I3, and of the hypercharge, Y (which is the sum of baryon number plus strangeness). Total isospin is 
approximately conserved when to„ = uid (and electromagnetic effect are ignored) and can be used to label baryon 
states in our world. However, when light quark masses differ arbitrarily, there is no conserved SU{2) subgroup of 
SU{3), and total isospin loses its utility. All of the baryons at the periphery of the SU{3) weight diagram in Fig.[7]have 
different I3 and Y and therefore do not mix with one another for any quark masses. However, the two states at the 
center of the octet weight diagram have the same I3 and Y, and can mix. In our world, the resulting eigenstates — the 
S° and the A — are approximate eigenstates of total isospin. In other worlds they will not be close to eigenstates of 
any SU{2) subgroup of SU{S) unless two light quarks have approximately equal masses. In general we label the mass 
eigenstates at the center of the weight diagram Chigh and Ciow 

It is a straightforward exercise in first order perturbation theory to compute the effect of the Hamiltonian, Eq. pil.3[) 
on the spectrum of octet baryons. For completeness, the analysis is presented in Appendix [A] and its results are 
summarized in Table [D The key point is that the Wigner reduced matrix elements (WRME's), conventionally labeled 
F and D, can be extracted from octet baryon mass differences in our world and provide the necessary dynamical 
information for computing octet baryon mass differences in any of the other worlds we study. 



B. Fitting parameters to the baryon mass spectrum 

In order to extract the values of the parameters F and D, we fit the matrix elements of iJftavor to the measured octet 
baryon masses after subtracting the electromagnetic mass contributions. In our universe a/ m§ + rri^ = 1.0002 mg « 
mg, so the quark masses 7713 and mg, always appear in combination with the WRME's F and D. We choose 

F 7778, Dms, — , and = (Fqcd) + m-o (Bq) (III.5) 

TTig 

as our fit parameters. 

There are eight baryon masses to be fit in terms of our four parameters. One linear combination of baryon masses 
is the Gell-Mann-Okubo relation and is independent of our parameters. Aq is simply the average octet baryon mass. 
This leaves six mass differences to be fit by three parameters, nisF, mg,D, and nis/ms- To keep a strong connection 
to chiral dynamics we take the ratio 7713 / TTig from the analysis of pseudoscalar meson masses rather than from our fit 
to baryon mass differences. According to the Particle Data Group (PDG), 

nil 777 

— = 0.56 ± 0.15; — = 20.1 ± 2.5 , (III.6) 

m-d ma 



The subscripts for the generators follow the usual convention for the Gell-Mann matrices. 
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Parameter 


Fitted value 




1150.8 MeV 


Fran 


-109.4 MeV 


Dms 


35.7 MeV 




0.0197 



TABLE II: SU{3) linear perturbation theory parameters. 



from which the PDG extracts mg/mg = 0.0197 ±0.0071 MeV [Tl|. If we fit mg/mg, we obtain mg/mg = 0.0181 MeV, 
which is quite close to the value quoted by the PDG. Note that we have not attempted to do a careful error analysis 
since the uncertainties in our theoretical modeling assumptions exceed the uncertainties in the parameters that we 
have extracted from data in our world. The values of the parameters of Eq. (jIII.5P that we use in our analysis are 
given in Table [III 



C. The matrix element of ?no(0o) and the pion-nucleon a term 

Analysis of flavor SU (3) symmetry breaking cannot give us information about cither the average quark mass or the 
matrix element of the SU{3) invariant operator 0j. However, SU{2) x SU{2) chiral symmetry violation can. It is well 
known that a measure of chiral symmetry violation in the baryon sector known as the pion-nucleon sigma term, CTttTV, 
can be extracted (with considerable theoretical input) from experiment, and can be related to the matrix element of 
mo0o- This is not the place for a critical review of the (still somewhat controversial) analysis of the chiral dynamics 
and pion nuclcon scattering that leads to a value for aj^N Instead we briefly quote the results of the standard 

analysis of Gasser and Leutwyler in [l3| and refer the reader to their review for further details. The extraction of the 
a term from low energy ttN scattering remains controversial. For more recent evaluations, see [13] and 

The aj^N term is defined as the matrix element of the SU{2) x SU{2) violating scalar density, 



{N\uu + dd\N) 

(iV| m„uu + rriddd \N) — —{niu ~ md) {N\ uu — dd \N) 
{N\muUU + mddd\N) 



(III.7) 



where |iV) is shorthand for the average of proton and neutron matrix elements. In going from the second to the third 
line, we have dropped a term that is second order in SU (2) violation, and therefore is negligible in our world. The 
combination m^uu-\-mddd can be written as a linear combination of the operators toqQo ^'^'^ Since the matrix 

elements of the latter are determined by baryon mass differences, knowledge of the cr-term fixes the matrix element 
of moOj. 

Rewriting Uj^n in terms of irreducible tensor operators, we have. 



V2mo + mg 
(^ttN = 7. ( N 



V2Ql + Q\ 



N 



(III.8 



Using the known matrix elements of 0g we obtain 



V2^ 

mg 



f ^(A^lmoeJiiv)- 

3 mo 



D 



(III.9) 



Given a value for (Ttttvi we can use the i^mg and Dm^ fit values and niQ/mg derived from the PDG quark mass 
ratios to solve for the singlet flavor mass term (A^| jtiqOJ thereby providing the needed separation of the average 
mass of the octet baryons, Aq, into a QCD component that does not scale with quark masses and this flavor singlet 
component that scales with the sum of quark masses rriQ. 

Gasser and Leutwyler conclude that the value of a-^N to be used in Eq. (jlll.Qp is aT^N = 35 ± 5 MeV. Substituting 
into Eq. (|III.9p and keeping track of the experimental uncertainties, we obtain 



(Ql) = 368 ± 101 MeV 



(III.IO) 
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FIG. 8: A slice through the quark mass prism at fixed xt = = 100. The axes X3 and xg, as well as xt are defined in Eq. IIII.13I Our 
universe corresponds to the point marked "us," near the lower right-hand corner. 



Using the values of the matrix elements of the other scalar densities, we can separate the contributions of the 
average u and d quark masses and the s quark mass, 

{N\m® {uu + dd)\N) = 35±5MeV 

{N\mfss\N) = 124±94MeV, (III.ll) 

which were the values quoted in Eq. pi.l3|) . 

Having determined mo(Oj), albeit with large errors, we can separate the QCD dynamical mass out of the the flavor 
independent term in (i?), 

(Hqcd) =Co^Ao- mo(ei) = 783 ± 101 McV (III.12) 

as quoted in Eq. (jII.13|) . Although the uncertainties in {N\moQQ\N) , and especially in {N\mfss\N) are large, the 
fact that (Hqcb) and (mo0Q) always appear in the combination Aq makes the analysis relatively insensitive to these 
uncertainties. 

It is well known that higher-order corrections in to^ are important to the baryon mass spectrum in our world. Such 
non-linearities mostly affect the flavor-singlet contribution to the baryon masses — as evidenced by the fact that the 
Gell-Mann-Okubo formula is accurate to within 1% — and are therefore primarily relevant to the issue of extracting 
of the (Ttt Af-term l43j . Since we have set up the problem for this investigation in such a way that we are largely 
insensitive to the uncertaintly in the value of cr^jv , nonlincaritics in mg will not affect our results significantly. Notice 
also that the central values for the quantities in Eqs. (jlll.llj) and ()III.12p agree with those given by Borosoy and 
Meissner [i^l, who take non-linearities in nis into account. 



D. Coordinates in the space of light quark masses 



Let us introduce dimensionless coordinates xt, X3 and 2:8 to describe the space of quark masses as pictured in Fig. 
2) To sidestep the problem that quark masses are not very accurately measured (and that they are rcnormalization- 
scale dependent), we will always work with quark masses and quark mass combinations as ratios to m® (the sum of 
the values of the light quark masses in our universe). With respect to the toqi 'tisj and mg defined in Eq. pil.4p : 

100(to„ + nid + rUs) 

XT = Tot X 



X3 



2m3 
^/3 



100 


= V6?T!,o X 


100 


100 


_ 100(m„ 


- rud 


V3m^ 


100 

m® 


100(to„ 


+ rud 




3m® 



_ 2ms 

" TT " = ^ • 

For a given xt {i.e., for a fixed sum of the light quark masses) the resulting two-dimensional space of quark masses 
may be pictured as the interior of an equilateral triangle with altitude xt. Choices of quark masses may be described 
by the orthogonal coordinates x^^g, as shown in Fig. [8l with ranges 



(III. 14) 
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Baryon 


Ct 




C3 


(-^em) 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


P 


3.68 


3.53 


1.24 


0.63 


n 


3.68 


3.53 


-1.24 


-0.13 




3.68 


-2.84 


2.44 


-0.07 




3.68 


-2.84 


-2.44 


0.79 


E+ 


3.68 


-0.69 


3.68 


0.70 


E" 


3.68 


-0.69 


-3.68 


0.87 


C*high 


3.68 


0.69^x1 + xl 


-0.21 


C'low 


3.68 


-0.69^x1 + xj 


0.21 



TABLE III: Values of the parameters in the equation baryon mass equation Mb = Co + ctXt + cgxs + 032:3 + {B\H-em\B) , for 
the quark mass coordinates xt, x^, and are defined in Eq. piI.13|l . As xt,3,s vary, the QCD contribution, Co, is adjusted to 
keep the mass of the hghtest baryon flavor muhiplet flxed at 940 MeV. Note that C§ = 783 ± 101 MeV. 



For our world, x® = 100, x® = -1.17 ±0.43, and x® = -59.5 ± 1.1. 

In order to express the baryon masses as functions of the scaled variables, xy, 0:3, and Xg, we must separate the 
flavor SU{?>) invariant contribution to the baryon masses, M{xt), into the piece independent of quark masses, Co, 
and the term proportional to {m^Ql), as defined in Eq. (|III.12|) : 

M{xt) = ^0 = Co + ^{N\m®Ql\N) = Co + ctXt (III.15) 

We can now express the mass of an octet baryon B as the sum of five terms: 

Mb = Co + ctxt + cgccg + C3.T3 + {B\H^m\B) . (III.16) 

The coefficients ct, C3, and cg and (H-em) are given in Table Hill Although both Co and ctXt have large uncertainties, 
the fact that Co (via its dependence on Aqcd) is adjusted throughout the parameter space to keep the average mass 
of the lightest flavor multiplet fixed removes the uncertainty from the parameterization of hadron masses. Therefore 
the errors for Co and (mo 6 J) are not quoted in the table. 



E. Baryon and meson spectra for a fixed sum of quark masses 



P=1/3 



P=2/3 




P=0; P=1 



As we discussed in Sec. IIIDl a triangular slice through 
the niu^d.s mass prism corresponding to a fixed value of xt 
will have a central region of type (i), where flavor SU{3) is 
a good symmetry and all of the octet baryons can appear in 
nuclei. If xt is large enough, then along the perimeter of the 
triangle (i.e., where the quark mass differences are greatest) 
there will be regions of type (ii) and (iii), where only the 
lightest pair of baryons can participate in nuclear physics. 

To understand the rescaling that defines our slice through 
the SM and some of the physical phenomena we have identi- 
fied, it is useful to examine the way that baryon masses vary 
over a triangular slice. First we track the baryon spectrum 
as we move along the perimeter of the triangle correspond- 
ing to XT = X® = 100. Figure [5] shows how we define the 
coordinate P to indicate where along the perimeter wc arc. 

Figures [TOl and [TT] show the octet baryon masses computed from Eq. (|III.16p (without the electromagnetic corrections), 
as functions of P. 

The baryon masses in Fig. [TO] are computed with a fixed value of Aqcd = ^qcd- However, as explained in Sec. 
II C[ we adjust Aqcd to keep the average mass of the lightest pair of baryons fixed at 940 MeV everywhere along 
our slice through the parameter space of the SM. Figure [11] shows the mass distribution after the rescaling of Aqcd ■ 
Candidates for congenial worlds occur when the two lightest baryons have nearly equal mass. As the figure shows, 
this occurs in two characteristic circumstances, first at P = 0, 1/3 and 2/3 (type (iii) regions in Fig. [5]), and second 
at P = 1/6, 1/2, and 5/6 (type (ii) regions in Fig. [5]). We will discuss these areas in detail in Section FVl 



FIG. 9: The definition of the perimeter coordinate used in 
Figs. [Tol through 1131 is shown graphically. Regardless of the 
value of xx, the coordinate P is defined to run from zero to 
unity around the triangle. 
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FIG. 10: Baryon masses (without electromagnetic corrections) at points along the perimeter of the quark mass triangle with = a;® = 
100 (i.e., for the sum of the three light quark masses equal to what it is in our world) and with fixed Aqoq = A^j^^. Sec Figurc[9]for a 
definition of the P axis. 



Mbaryon(MeV) 

A 

1,400- 




FIG. 11: The same baryon masses plotted in Fig. 1101 are shown again, but now adjusted by an appropriate rescaling of Aqcd so that 
the average mass of the lightest pair of baryons is held fixed at 940 MeV. For the baryon masses shown in this plot, Aqcd varies between 
1.03A®„„ and 0.85A®„„. 



The formulas for the baryon masses are Unear in and xg (or in the radius \/ x'^ + a;|), so considering a smaller 
triangle would yield similar but rescaled plots, with the curves in Figs. ITOl and fTTl squeezed together. This is true of 
both the perimeter of a triangle that corresponds to a smaller value of xy, as well as of a triangular trajectory interior 
to the original triangle. Two such triangles with the same baryon spectrum are shown in red in Fig. 1121 

The adjustment of Aqcd to keep the average mass of the two lightest baryons fixed at 940 McV eliminates the 
dependence of the baryon spectrum on xt, but even though the two triangles in red in Fig. [T^lgive the same baryon 
masses as one goes around their perimeter, their physics is not equivalent because the meson spectrum is not the 
same in the two cases. To see this, consider the masses of the light (flavor SU (3) octet) pseudoscalar mesons. From 
chiral perturbation theory: 

ml Bo{mu + nid); to^± = i?o("i„ + m^); to^o = -Bo(TOd + m^) , (III.17) 

where Bq is proportional to Aqcd- These vary as functions of xt.3.8 because of the change in the quark masses and 
also because of the change in Aqcd required to keep the lightest pair of octet baryons at an average mass of 940 
MeV. Notice, though, that the two central states of the meson octet, tt" and rj, mix. We label the corresponding mass 
eigenstates as ^high and ^low, with masses: 



™lhigh/iow ^ ( ± Jm1 + m'^ + m'^ ~ ruumd - mums - mams j . (III. 18) 
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Figure [13] (a) shows the spectra of baryons and mesons, as 
a function of < P < 1/6, for the perimeter of the xt = 100 
triangle. The other sectors in P are obtained by symmetries. 
For clarity only the three lightest baryons are plotted. Fig- 
ures [13] (a) and (b) show, respectively, these spectra for the 
perimeter of the xt = 40 triangle, and for a trajectory inside 
the XT = 100 triangle that corresponds to the projection of 
the XT = 40 triangle. The baryon spectra are identical in 
these last two cases (because of the way we have defined our 
slice through the parameter space of the SM), but the meson 
spectra are somewhat different. As previously mentioned nu- 
clear binding may depend more strongly on the mass of the 
/o(600), which is likely to be less affected by changing quark 
masses (especially when Aqcd is readjusted as we have de- 
scribed), than on the masses of pseudoscalar bosons. So we 

arc optimistic that our analysis can be applied over a relatively wide range in xt as long as low energy nuclear physics 
docs not change qualitatively, as it does when more than two baryons appear in stable nuclei. 

In Fig. 1131 a vertical line identifies the value of P for which the mass of the pion is the same as in our world. This 
corresponds to the position where the m^r = 140 MeV line shown in Fig. [HUa) cuts the perimeter of the corresponding 
triangle. Around that value of P, the nuclear interaction between the two lightest baryons (the proton and the neutron) 
would be very similar to what it is in our world. For P = 1/6 the meson spectrum is qualitatively quite different. 
The lightest meson is a mixture of an isosinglet (ry) and the /a = component of an isovector (tt*^). Furthermore, the 
neutral kaons and the charged pions have become degenerate in mass and are only slightly heavier than the lightest 
meson. It would be very difficult to describe precisely the features of the nuclear interaction between the lightest 
baryons (the neutron and the S") in such a world. 



FIG. 12: The perimeter of a triangle of small constant xt is 
shown projected onto a central region within a triangle corre- 
sponding to a larger value of xj'. The two triangular trajec- 
tories, drawn in red, give the same spectrum of octet baryon 
masses, but different meson spectra. 
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FIG. 13: The masses of the three lightest octet baryons and of the the light pseudoscalar mesons are shown as functions of the perimeter 
coordinate P, for < P < 1/6. The vertical line identifies the point where the pion mass is the same as in our world. The spectra are 
shown for: (a) the perimeter of the xt = 100 triangle, (b) the perimeter of the xt = 40 triangle, and (c) a trajectory within the xt = 100 
triangle that corresponds to the projection of the perimeter of the xj- = 40 triangle. 
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IV. MODELING NUCLEAR INTERACTIONS IN OTHER WORLDS 

In order to explore the congeniality of worlds with different quark masses we need to at best compute or at least 
estimate the binding energy of nuclei in those alternative worlds. In worlds where nuclear forces are the same as ours 
but the masses and charges of participating nucleons are different, we can compute nuclear binding energies of low 
atomic mass nuclei accurately by comparing with analog nuclei in our world [19[. This method enables us to judge 
the congeniality of all the worlds where two light quarks participate in building nuclei. To study higher atomic mass 
nuclei in these worlds and to estimate the properties of nuclei in worlds with more than two participating baryons, we 
turn to a generalization of the Weizsacker semi-empirical mass formula (SEMF). As noted in Sec. II Al we have only 
succeeded in making qualitative statements about worlds with three very light quarks. Nevertheless the generalization 
of the SEMF to these worlds is interesting in itself and appears to be new, so we include it here. 

A. Analog nuclei 

As explained in Section II, we can explore a wide class of interesting worlds while keeping the dynamics that 
determine nuclear interactions approximately fixed. Nuclei in these worlds differ from ours in two important ways: 
First, the charges of the two participating baryons may differ from ours. If the s and d quarks are lightest then the 
5]~ and 2~ replace the proton and neutron. If the c and u are lightest, then the S^"*" (uuc) and the (ecu) are the 
constituents of nuclei. Second, the masses may differ. However, the average mass of the nuclear constituents has been 
kept at 940 McV and the mass difference ranges only over ±20 — 30 MeV. Altering the charges affects only the nuclear 
Coulomb energy, a change that is easily made in the framework of the semi-classical estimate of the Coulomb energy, 
= (0.7MeV)Z(Z- l)/Ai/3. Altering the masses does not affect the nuclear binding energy. It does, however, 
affect nuclear decay patterns since it determines the thresholds for /?-decays and for weak nuclcon emission. 

It remains only to determine the appropriate analog nucleus for the case of interest and make the necessary mass 
and charge adjustments. As a simple (but important) example, consider the stability of "oxygen" in a world where 
the s and d quarks are lightest. In that world the nucleus with the chemistry of oxygen has the baryon content 
(S~)^(S~)*. It is helpful to introduce a new notation to describe the nuclear and chemical properties of atoms in 
worlds where the constituents of nuclei have different charges than in our world. If we limit ourselves to worlds with 
two participating baryons, then an atom can be labeled by the number of nucleons of each type, A^i and N2, the baryon 
number, A ~ Ni + N2, and the electric charge, Z , which we denote by the chemical element symbol independent 
of whether the charge is positive or negative {e.g. Z = ±4 =J> Be). For convenience, we will represent the form of 
element El with A^i baryons of type 1 and N2 baryons of type 2 by the symbol -^_^F,\j^^. Thus chemical oxygen in a 
d-s world would be IO4. This oxygen has the nuclear binding of 4Be4 in our world, plus some additional Coulomb 
energy. Since 4Be4 is unbound in our world, and the additional Coulomb energy in a d-s world only unbinds it further 
relative to its decay products (2262, the (S~, S~) analog of the tightly bound a particle), we conclude that chemical 
oxygen is not stable in s-d worlds, regardless of the s-d mass difference. 

In worlds that differ from ours only in the values of the u-d quark mass difference, the analog method works 
somewhat differently. The binding energies of nuclei are the same as in our world, however the systematics of weak 
decays change. It is straightforward to determine the limits of nuclear stability. The deutcron provides a simple 
example. It is bound by B(^H]^) = 2.22 MeV in our world. Changing the mass of the proton and neutron does not 
alter this. However when Mp — Mn > Hj^) — rUe = 1.71 MeV, the deuteron can decay to two neutrons by electron 
capture (^He~ nnVf.). Alternatively, when M„ — Mp > B(^H]^) + jTig = 2.73 MeV, the deuteron /3-decays to two 
protons (^H— > ppe~9e)- These are the simplest examples of weak nucleon emission, which figures importantly in 
determining the limits of congenial worlds. 

B. Generalizing the Weizsacker semi-empirical mass formula 

Weizsacker's semi-empirical mass formula (SEMF) provides an excellent effective Hamiltonian for the nuclear ground 
state, when the number of each participating baryon species is conserved. It is usually written as a formula for the 
nuclear binding energy, but with a change in sign and addition of the nucleon rest masses, it describes the total energy, 
which is better suited to our purposes. In our world, 

Hm}) = ^^M^ ~ ^ + ^f^y^ + ^f^^^ + + (IV.l) 
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Parameter 


Value (MeV) 


© 


15.56 


e 


17.23 


e® 


0.7 


e® 


93.12 


A ffi 

A* 


12 


^© = ^r© _i_ 

Cv — '^v ^ 22/5 
o8/5 e 

— "^a 


32.0 
39.8 
60.8 



TABLE IV: Parameters in the SEMF, as quoted in [31 . The first five are the standard parameters as defined in Eq. (|IV.l|l . 
The last three are the parameters that enter the Fermi-gas modified SEMF developed in Sec. lIV]T2l Note that the superscript 
® is a reminder that these parameters are fit to nuclear masses in our world. 



where A = ^ A^,; = Z + A^, /a = (Z — N)/2, and as before, the superscript © denotes parameter values in our world. 
The volume, surface, Coulomb, and asymmetry terms are self-explanatory. The last term is the pairing energy, j] is 
+ for odd-odd, — for even-even, and for odd-even and even-odd nuclei. We take the values of the coefficients, e®, 
from [Hi. They are listed in Table [IVl 

To consider different nuclear constituents, we must generalize the SEMF in two ways. First we must allow for 
different symmetry groups. As we shall see the SEMF assumes isospin symmetry, and if we are to consider worlds 
with three light quarks we must generalize to SU(3) flavor symmetry. Second, although isospin or SU{3) may be an 
excellent symmetry of the strong interactions if quark masses are small enough, the symmetry may be dramatically 
broken by kinematic effects when a baryon becomes so massive as to drop out of nuclear dynamics. This effect requires 
some explanation. Consider a sequence of universes where the u and d quark masses are fixed at the values in our 
world, and the s quark mass is slowly increased from a few MeV. For the smallest s quark mass, all eight octet 
baryons would be expected to participate in nuclear physics and SU (3)-flavor symmetry rules. When the s quark 
mass increases to the point that the lightest strange baryon mass exceeds the proton and neutron Fermi energy in 
the nucleus, strange baryons disappear from nuclear dynamics. This occurs when the hyperon masses are only a few 
tens of MeV heavier than the nucleon mass, in a regime where the SU (3) symmetry of the strong interactions is still 
quite a good symmetry. We believe that this essentially kinematic, but massive, violation of flavor symmetry can be 
incorporated into the SEMF by modeling part of the volume and asymmetry energies by a degenerate Fermi gas. 



Restoring flavor symmetry to the SEMF 



Imagine that the proton and neutron masses were exactly the same, so the strong interactions were exactly SU (2) 
invariant. Two terms in our world's SEMF, the Coulomb and asymmetry terms, appear to violate SU{2) flavor 
symmetry. The Coulomb term violates SU (2) because it originates in electromagnetic interactions which, of course, 
are not SU{2) invariant. However, what are we to make of the appearance of the explicitly SU(2) symmetry violating 
operator oc /| in the asymmetry term of Eq. pV.l[) ? The asymmetry energy is a strong interaction effect and must 
depend only on SU{2) flavor invariant operators. 

The solution to this puzzle is that /| should be replaced by the quadratic Casimir operator of SU{2), I{I + 1), and 
that electromagnetic violations of SU{2) allow I to adjust to its minimum value, which is I = I3. In the exact SU{2) 
limit and without electromagnctism, nuclei with a given A would form degenerate multiplets labeled by the total 
isospin I. The lightest multiplet would have the lowest value of /. Since I > I3 and I3 is conserved in the absence of 
the weak interactions, this occurs when I = I3. 

Thus, electromagnctism allows the states with / > /a to decay quickly to the ground state, so for all intents and 
purposes the S'C/(2)-symmetric asymmetry term, /(/ -I- 1)1/=/^^;,^ can be replaced by /3(/3 -I- 1). Switching between 
this and the usual form, /|, produces a negligible shift in the computed values of the binding energy B{A) and the 
charge Z{A) of the lightest nucleus with given A. That is, the two forms of the SEMF have nearly the same valley of 
stability. Therefore, a more illuminating way to write the SEMF close to the S'C/(2)-flavor symmetry limit, is 



A'^''^ A 



'^1/2 



(IV.2) 
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where the superscript on reminds us that this effective Hamiltonian apphes near the SU (2) hmit (regardless of 
which two quark species are Hght). 

The same analysis must be performed near the SU{3) symmetry limit. The parameters in the SEMF are no longer 
the same as in our world (so we replace e® and A® by unknown ci and A), the variable r] must be redefined to count 
the number of unpaired baryons, and the operator that appears in the asymmetry term must be the quadratic Casimir 
operator of SU{3), which we denote by G^. The value of depends on the representation of SU{3) to which the 
nucleus belongs. SU{3) representations can be labeled by two parameters, X = ni — n2 and fi = n2 — n^, where 
ni, 712,713 are the number of boxes in the rows of the Young Tableau that defines the representation. The expression 
for G2(A,Ai) is 

G^iX,fi)^\{X'' + Xfi + fi^)+X + li (IV.3) 

and the S't/ (3)-SEMF can be written as 

Z{Z-1) G^{X,n) 



.41/3 ^ A 



(IV.4) 



In the S't/(3) case, the electromagnetic and strong interactions will quickly adjust the flavor content of the nucleus 
to minimize the asymmetry energy. In contrast to the SU{2) case, strong interactions play an essential role here. 
In the two flavor case the two conserved quantities, A and I3, uniquely fix the number of both types of baryons, 
which therefore remain fixed as / adjusts to its minimum. With three fiavors there are as many as eight participating 
baryons, but only three conserved quantities, A, I^, and the hypercharge, Y. So the strong and electromagnetic 
interactions together will adjust the abundance of the eight species to minimize the ground state energy with fixed 
A, h and yis. 

We need, therefore, to compute the minimum value of G^(A, /i) at fixed values of ly, and Y . Fixing /a and Y fixes 
a point in the SU{Z) weight diagram (labeled as usual with along the horizontal and Y along the vertical axis). 
A representation is specified by a hexagon (or when A or /.i vanishes, a triangle). The values of A and /i give the 
lengths in J3 of the topmost and bottommost horizontal lines in the diagram. It is easy to see that the value of 
at fixed I3, and Y is minimized when the point (/a, Y) appears on the boundary and at a corner of the polygon. The 
G^(Aniin7 Mmin) resulting from this procedure can be explicitly computed: 

G2(A,/^)|._ ^ll+'^Y^ + \ {\Nu - Ni\ + \Nd-Ns\ + \N, - iV„|) . (IV.5) 

Here Nu,d,s are the coordinates of the point in question along axes that count the number of u, d, and s quarks. These 
terms can be written as piecewise continuous functions of and Y . Whether or not they can be ignored depends on 
how close the system is to the symmetry point /a = F = 0. In Sec. IVE 21 where we use this formula, we keep these 
terms. Although they make a significant quantitative difference, they do not affect the physical conclusions in that 
section. Eq. (jIV.4p and Eq. (jIV.SP together provide the SU (3) flavor generalization of the SEMF that we seek. 



2. Kinematic flavor symmetry violation in a Fermi gas model 



All of the qualitative features of the SEMF except for the pairing energy can be accounted for by treating the nucleus 
as a non- interacting Fermi gas of baryons under constant pressure. Furthermore the Fermi gas model includes, in a 
physically intuitive way, the symmetry breaking effects that are essential for our analysis of worlds with different quark 
masses. In particular, if the mass of a given baryon species exceeds the Fermi energy, then the species disappears from 
the nucleus. This wholesale flavor symmetry violation must be taken into account even when the dynamical breaking 
of flavor symmetries is small. 

The non-interacting, degenerate Fermi gas under constant pressure gives a one parameter model for nuclear energies. 
If the pressure is chosen to reproduce the measured Fermi energy in heavy nuclei, then the surface. Coulomb and 
asymmetry energies are predicted. The surface energy arises from the depletion of the density of states near the 
surface of the nucleus. The Coulomb energy is that of a uniformly charged sphere whose volume is fixed by the 



The reader will notice that the appearance of the baryon numbers {Ni} in Eq. |)IV.4|I is problematic since they are dynamically 
determined (though constrained by the values of A, I3 and Y). Below we shall explain how they are fixed in a Fermi gas model. 
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pressure and the atomic mass. The asymmetry energy arises because the kinetic energy is at a minimum when the 
proton and neutron Fermi momenta are equal and grows when they differ. However, the degenerate Fermi gas model 
accounts for only a fraction of the actual effects as parameterized in the SEMF. This is not surprising: some of the 
nuclear surface energy comes from the different interactions of nucleons on the surface relative to the interior of the 
nucleus, and some of the asymmetry energy comes from the dependence of the nuclear forces on isospin (the simplest 
example being the difference in the / = (deuteron) and 1=1 channels of proton-neutron scattering). 

In order to model flavor symmetry violation in worlds where two baryons participate in nuclear physics we construct 
a degenerate Fermi gas model with the following ingredients. (We address the three-flavor case later.) 

• We seek an effective Hamiltonian as a function of the masses and numbers of each of the baryon species. When 
we wish, we can flnd the most stable species at each A by minimizing over Ni at fixed A. 

• The kinetic energies of the baryons are given by a degenerate Fermi gas under constant pressure, with the 
pressure chosen to reproduce the Fermi momentum in heavy nuclei in our world. 

• We take the surface and Coulomb energies directly from the SEMF, ignoring possible flavor symmetry violating 
effects in these terms because they are higher order. 

• We compute the Fermi gas contribution to the asymmetry energy and subtract it from the phenomenological 
asymmetry energy term in the SEMF. The remainder we ascribe to interactions, which we assume to be flavor 
independent because we work in worlds where SU (2) symmetry violation is small. 

• The important, dramatic flavor symmetry violation is thus contained in the Fermi gas kinetic energy. 

Let Ni be the number of baryons of the i'^ species in the nucleus and ki be its Fermi momentum. Ni and ki are 
related by 



V 



(IV.6) 



where the volume V is yet to be determined. The total internal energy including the baryons' rest mass contribution, 
is 

i i ^ ^ 

The i^^ chemical potential is 



Mi 

and the pressure is 

P = 



dU 

mi 



dU 

av 



V 



Note that we consistently keep baryon mass differences in the rest energy, where they drive large symmetry violating 
effects, but ignore them in the kinetic energy, since they enter as AMi/M"^ ~ 10~^/M, which is negligible. We use 
M to stand for the average mass of the (two) participating baryons. 

It is convenient to use the grand canonical ensemble to describe a gas under constant pressure. The grand potential 
is defined as a function of {iii} and V , 

i i 

and its intensive counterpart is 

c..(^,) ^ ^M^flll) ^ __J_ [2M{^,. - M.)f' . (IV.ll) 



The quantity Q is constructed so that 

dfl 
dfii 



= -N, . (IV. 12) 
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The total energy is the internal energy of the Fermi gas together with the work it does against the confining pressure 



Dynamical equilibrium is achieved by minimizing E with respect to V at fixed Ni, 



dE 

w 



E 



d^i dV 



Ni 



dfM 
dV 



Ni 



Po ^ . 



But from Eq. pV.12p . 



V 



^dujj 



dfii 



-N,, 



and, substituting back into Eq. ()IV.13P 



^uj,{fi,) + Po = , 



E{{N,})^J2^'^^^ 



(IV. 13) 

(IV. 14) 

(IV.15) 
(IV.16) 

(IV. 17) 



Equations (|IV.16p and pV.17p are the fundamental dynamical principles for the system. Using them it is possible 
to solve for all of the dynamical properties of the system (energy, volume, compressibility, etc.) in terms of Pq and 
the baryon numbers, {Ni}. Of course, we are interested in the total energy, which is 



.5/3 



3/5 



with 



1 



2A/3/5 



(IStt^Po)'/'. 



(IV. 18) 



(IV.19) 



We fix the value of Po based on the Fermi momentum of nucleons in heavy nuclei in our universe. Quasi-elastic 
electron scattering on heavy nuclei yield values of fcp = fc^ = 245 MeV [43], which yield a pressure P® ~ 0.827 
MeV/fm3 and e[f = 32.0 MeV. 



3. Modified SEMF in the two flavor case 

In this section we consider the two flavor case with two participating baryons with A = and /a = {N1—N2) /I. 

Equation (|IV.18p contains the kinetic contribution to the asymmetry energy which we seek. In the usual SEMF, 
however, this contribution is inseparable from the dynamic contribution which has been fit to the form ea^f/^ for 
small asymmetries. To determine the dynamical contribution, we must expand E{Ni, N2) around A'^i = iV2, identify 
the term with the form of the asymmetry energy, and subtract it from the phenomenological fit, leaving the dynamical 
term in the asymmetry. Expanding E a.s a. power series in I3, we obtain 

E{Ni,N2) = NiMi + N2M2 + ^^oA + ^eo^ + 0{ll/A^) (IV.20) 

When we explicitly include the Fermi gas kinetic energy we must compensate by removing the term proportional to 
A in Eq. (|IV.20p from the volume energy and the term proportional to « /3(/3 + 1)/^ from the asymmetry 

energy. Putting this together we obtain an SEMF 



3/5 



Z(Z-l) A 
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which incorporates the flavor symmetry violating effects that come from the Fermi gas kinetic energy. As long as the 
low-energy nuclear interactions are similar enough to those in our world, the e^'s in Eq. pV.2ip can be replaced with 
ef's and A can be replaced with A®. 

It is convenient to define modified coefficients of the volume and asymmetry terms, 

e 

= 4- ■ 

V W 22/5 ' 

-e^ EE ^ . (IV.22) 

These parameters — along with e® — are listed in Table IIVI 

In the context of the Fermi gas model, the barred parameters represent the dynamical contributions to the volume 
and asymmetry energy in contrast to the kinematic contributions that are contained in the term proportional to e®. 
In terms of these parameters the modified SEMF becomes, 

- Mh + l) j, Z{Z-l) A^ 



^. Modified SEMF in the three flavor case 

The analysis of the previous section can be generalized to the three flavor case by replacing the operator in the 
asymmetry term by the operator we derived in Section IIVB 1[ 

4=1 4=1 

where, for simplicity we have dropped the small corrections to the quadratic Casimir operator that are linear in the 
quark number differences. These terms can be included with little increase in complexity. Note, also, that we have 
removed the superscript "0" on the parameters as a reminder that we do not know the values of the SEMF parameters 
in the SU{3) limit because the physics of nuclear binding is qualitatively different than in the two flavor case. For the 
same reason, there is no point in trying to separate Ca h^to the dynamical and Fermi-gas kinetic pieces since neither 
is separately known. 

Eq. (|IV.24|) neglects an important piece of the dynamics that appears when we try to generalize the SEMF from 
two to three flavors. In the two flavor case the two quantum numbers conserved by the strong and electromagnetic 
interactions, A and I3, are sufficient to fix the numbers of both nuclear species. In the three flavor case, fixing A, I3, 
and Y does not fix the number of all eight baryon species. Strong interaction processes like 

S" + A ^ 39 + S" (IV.25) 

will establish "chemical" equilibrium among the various baryon species in order to minimize the energy of the system. 
This is equivalent to minimizing H^^^ {{Ni}) with respect to the Ni at fixed A, I3, and Y. In general this minimum 
cannot be found analytically even in the simple Fermi gas model. However in Section |V] we will argue that the case 
when the baryon mass differences are small compared to the scales of nuclear physics (eo, Cv, etc.) is of particular 
interest. In this case the minimization problem can be solved analytically to first order in the mass differences, giving 
an effective Hamiltonian (still ignoring the linear terms in the SU{3) Casimir operator). 



H^'\A,h,Y) = H^'H{N,}) 



Ni=N°{Aj3,Y) 



MA + -^A + ^'l^'^}^ + \hAMj + \YAMy - e^A + e.A^/a 



\Z-1) A 

^1/3 +^'^1/2 



4^ + .^+,^. (IV.2«) 
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Here M is the average (over the octet) baryon mass, and AM/ and AAfy arc the isospin and hypcrcharge weighted 
baryon mass differences, 

8 8 

AA// = ^ M,h, and AA/y = ^ M^F, . (IV.27) 

i=l i=l 

As expected, if the baryon mass differences vanish, the minimum in H^^\A,h,Y) is at /a = y = and the 
expansion of the Fermi gas kinetic energy to second order in and Y yields terms of the same form as the volume 
and asymmetry terms. Collecting terms with the same functional dependence, effectively reversing the logic after 
Eq. (jIV.2f [) . we can write the SEMF expanded to leading order in baryon mass differences about the SUli) limit, 

if (3) {^A, hX) = MA + i/gAAf/ + iyAA/y - e^A + e^A^I^ 

+ ^a^-^ + ^c^^ + r;-^ , (IV.28) 

where and are defined in terms of ev, Ca, and eo by relations of the form of Eq. pV.22p . Equation (|IV.28p will 
allow us to explore the charge versus mass relationship of nuclei in worlds with very small quark mass differences. 



5. Nuclear stability 

For completeness we review the conditions for the various decay mechanisms that figure into our studies of conge- 
niality. In our world weak interactions — various forms of /3-decay — act to minimize the energy at fixed baryon number 
and strong interactions, notably a emission, lead to baryon number changing decays. In the worlds we consider these 
processes are augmented by weak baryon emission. 

The criteria for fixed baryon number weak decays are elementary. For decay, Mi > Mf + me, and for electron 
capture Mi + rrie > Mf, where i is the initial state and / the final one. If nuclei are made of baryons of equal 
charge, e.g. the S~) world mentioned in Sec. IIV Ai weak decays are somewhat suppressed because they would 
correspond to strangeness changing neutral currents. This suppression will be discussed further in Sec IV Bl when it 
arises. 

The criterion for strong particle emission {e.g. a decay) is that the final state be more bound than the initial 
state. We need only consider the most stable nucleus, denoted [A\, for a given value of A since we are only interested 
in absolute stability. Let 6max(^) be the binding energy per baryon of this nucleus, and let [a] denote the nuclear 
fragment {e.g. the a particle) that it might emit and h{a) its binding energy per baryon. When the baryon number 
is large enough that the binding energy can be regarded as a smooth function of A, the decay [A\ [A — a] + [a] is 
energetically allowed if 

W(A)+A ^^— <6(a). (IV.29) 

The other decay process important for determining boundaries of stability is weak nucleon emission, which does 
not occur naturally in our world, but becomes important at large baryon mass differences. In weak nucleon emission, 
a heavier baryon weakly decays into a lighter baryon while it is emitted (with an accompanying lepton pair as 
appropriate). The nucleus gains the energy from the mass difference between the baryons. Exactly such a process is 
responsible for the semi-leptonic decays of hypernuclei in our universe. 

Next, we compute the threshold for weak nucleon emission. For concreteness we imagine that the relevant baryons 
are the proton and neutron. There are two possibilities: First, if Mp > Mn then the proton in the nucleus can capture 
an atomic electron to give a neutron, which is emitted from the nucleus, and a neutrino, 

WEAK NEUTRON EMISSION: [Z] + ^ '"^""^1 [Z - 1] + ?1 + I^e , 

[-41 [Z] ^ [^-11 [Z -l]+n + e+ + 17^ . (IV.30) 

Positron emission is also possible, as described in the second line of Eq. (jIV.30p . but the threshold is higher, so we 
limit our attention to electron capture. Second, if M„ > Mp then the neutron in the nucleus decays to pe^v^ and the 
proton is ejected, 

WEAK PROTON emission: [Z] ^^^'^^[Z] + p + + ■ (IV.31) 
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For definiteness we consider the first case, weak neutron emission, and quote the result for weak proton emission at 
the end. Weak neutron emission is energetically allowed when 



M{N, Z)+me> M{N, Z - 1) + M„ 



(IV.32) 



Let B{A, I3) be the nuclear binding energy in terms of A = N + Z and I3 = \{Z - N). Let AAf = Af„ - Mp 
Eq. (jIV.32[) can be rewritten in terms of binding energies, 



B{A,h)-B{A-l,h--)<m,-^M 



(IV.33) 



For very light nuclei this condition would have to be considered on a case by case basis. For larger A we can 
approximate the discrete difference as a derivative. 



dB 
1)1 



1 dB 

2 a?3 



(IV. 34) 



We are interested in nuclei already at the minimum of the valley of stability set by the weak interactions, so I3 = -^3(^) 
is set by 



dM{A,h) 



dh 







(IV.35) 



Since M{A,Ij,) — ^A{Mn + Mp) — I^AM — B{A,l3), a nucleus at the weak interaction valley's minimum satisfies 



dB 



-AM 



(IV.36) 



Substituting this back into Ec^. (|IV.34|) . we find that the nucleus at the bottom of the valley of stability is unstable 
to weak neutron emission when 



dB 
dA 



< me 



1 



-AM . 



(IV.37) 



Finally we relate dB/dA\j^ to the derivative of the binding energy along the valley of stability, which is defined by 

B^,^{A)=BiA,h{A)) (IV.38) 
where hiA) is determined by Eq. pV.35[) . Differentiating i?max(^), we obtain, 

dhiA) 



dB„,,^{A) ^ dB 
dA dA 
dB 
dA 



I3 



dB 
dh 

AM 



A dA 
dh{A) 
dA 



(IV.39) 



from which we can substitute into Eq. (jIV.37[) , to obtain a condition for instability in terms of the the familiar function 



WEAK NEUTRON EMISSION: — - < 771^ ~ AM 

dA 



1 dhjA) 

2 dA 



(IV.40) 



The nuclear asymmetry energy keeps dh/dA small throughout the regime of stable nuclei. In our universe it is zero 
for nuclei at the bottom of the valley of stability up to yl ~ 40, and increases only slowly thereafter. It is of order 0.1 
for heavy nuclei. Therefore to a good approximation (ignoring also the small mass of the electron) the criterion for 
instability by weak neutron emission reduces to 



WEAK NEUTRON EMISSION: 



dBmax(A) . 1 



dA 



< --AM 
2 



(IV.41) 
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The instability for weak proton emission, which occurs when the neutron is considerably heavier than the proton, 
are 

WEAK PROTON EMISSION: , . < -TO„ + AM ' 



dA \2 dA 

< 1am . (IV.42) 

dA ^ 2 ^ ' 

In worlds with two light quarks and two baryons participating in nuclear physics, the binding energy per nucleon in 
medium to heavy nuclei is approximately 8 MeV. Therefore weak nucleon emission instability sets in throughout a 
broad region of A when |AA/| > 16 MeV. 

Our result in Eq. pV.42p differs by approximately a factor of two from the result claimed in Sec. 4.3 of 
[ist . There the authors compute the constraint on the proton-neutron mass difference from the condition that 
complex nuclei be stable against /3-decay. Their Eq. (27), which reads m„ — — TOg < £^bin, should read 
m„ — TOp — TOg < — clB/dlslj^, since the process they are considering takes place at fixed A. In any event, this 
is not a condition for the stability of complex nuclei, but rather the condition that sets the isospin (or, equivalently, 
the N/Z ratio) at the bottom of the valley of stability, as a function of AAI. The process that limits the stability of 
complex nuclei is weak nucleon emission, and it is a coincidence that our own result differs from the result of [48| by 
a factor of ~ 2. 



V. CONGENIAL AND UNCONGENIAL WORLDS 



Using the tools developed in Sees. IIIII and IIV| we now proceed to investigate which values of the quark masses 
(along our choice of a slice through the parameter space of the SM as defined in Sec. [TTI would correspond to congenial 
worlds. We shall consider in turn the case of 

(a) one light quark leading to a single light baryon, 

(b) two light quarks with equal electric charge, 

(c) two light quarks with different charge, 

(d) one light quark leading to two light baryons, 

(e) three light quarks. 

A complete characterization of (e) is left for future investigation. 



A. One light quark leading to a single light baryon 

A universe where only one baryon is light enough to have a chance of participating in nuclear physics would be 
uncongenial. For instance, if the d and s were so much heavier than the u that the decuplet uuu state (i.e., the 
A++) were lighter than all the octet states, then no complex nuclei would exist [T^]. A possible bound state of two 
A"'""'''s would have as its analog the dineutron, which is unbound due to its high asymmetry energy penalty. The 
electromagnetic repulsion would further obstruct binding between A^+'s. The only chemical element would be helium 
(if the light quark has charge 2/3) or hydrogen (if it has charge —1/3). 

For the lightest decuplet baryon to be lighter than all the octet baryons somewhere in the space of quark masses 
corresponding to a fixed my = + rrid + Wg, the value of my would have to exceed some lower bound, which can 
be estimated as follows: Consider the extreme case m„ = 0, ^ = mT/2, where the possibility of a light A++ is 
greatest. The lightest octet baryons would be the proton and the S+, both with mass fixed to 940 MeV on our slice 
through SM parameter space. From Eq. (jIII.16p and Table [iTT] (ignoring the electromagnetic correction) we then have: 

Mp = 940 MeV = (783 ± 101 MeV) + (2.71 MeV) xt , (V.l) 

^QCD 

where the factor 2.71 MeV comes from evaluating the coefficients in Eq. piI.16p when xs = —xt/6 and .T3 = — xt/2\/3, 
which follow when rriu = and nis = rud- 

To accurately estimate the mass of the A^"*" in a world with 7n„ = and rn^ ~ rris, we would need the decuplet 
analog of the cr-term (49j . For our purposes, however, it will suffice to make the quark model motivated assumption 
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that the mass of the A++ is independent of rud and and that its dependence on m„ is the same as the proton's 
(which, incidentaUy, agrees with the estimate of [i^). This gives 

777® 

(A++|777„?2w|A++) ~ ^ X 35±5McV w 25±4MeV , (V.2) 

so we estimate the mass of the A++ to be « 1210 MeV in a world with Aqcd equal to its value in our world, but 
with the u quark mass set to zero. In other words. 

Ma = 1,210 MeV ^§2£ . (V.3) 

^QCD 

Then the maximum value of xt such that the A^"+ remains heavier than the proton for ?ti„ ~ is 

Let us summarize the situation at fixed xt'- For worlds with xt well below a;™'''^, octet baryons are the lightest 
stable baryons throughout the mass triangle, and decuplct states can be ignored. As xt grows beyond x™'^^ = 122 ±29 
MeV, a decuplet state is less massive than the lightest octet baryon state in three regions at the perimeter of the 
triangle near the midpoint of each side (P 1/6, 1/2, and 5/6). These regions are uncongenial. This is an important 
effect because these regions (type (ii) in the classification of Fig. [5]) may otherwise have been congenial. We return to 
this point in Subsection IV Dl below. 



B. Two light quarks of the same charge 

Consider first the case of two light quarks with charge —1/3. We shall sec that these worlds arc uncongenial because 
it is impossible to make stable nuclei with charge greater than 4. To establish this it is sufficient to consider the case in 
which the light quarks are degenerate in mass. Otherwise, the nuclei at the bottom of the valley of stability (defined 
by the condition dM/dN2\A=Ni+N2 = 0, where Ni^2 are the numbers of each of the participating baryons) would be 
less bound than in the degenerate case, due to the asymmetry energy penalty. 

In these worlds the participating baryons arc the S~ and the S~. Because both baryons have charge —1, the charge 
of a nucleus. Z = —A. With no mass or charge difference between the baryons, the valley of stability lies on the line 

Our universe does have a small M„ — Mp quark mass difference, but it is much smaller than the Fermi energy and 
thus is negligible in determining the valley of stability. Hence in our universe (before the Coulomb energy becomes 
significant) we can find analog nuclei that lie on the corresponding line Z = A/2. 

Returning to the world with two light -1/3 charge quarks, we begin by considering very light nuclei for which the 
SEMF cannot be used. Using the notation introduced in Sec. IIV Al we denote the first nucleus, with A = 2 and 
Z = —2, as iHC]^. Its analog in our world is deuterium. We may estimate its Coulomb energy by modeling the 
deuteron nuclear potential as a finite spherical well with depth Vq and radius a. To determine Vq and a we impose 
the condition that the nuclear binding energy match its measured value for the deuteron, B = 2.22 MeV, and the 
r.m.s. radius correspond to the observed deuteron diameter, 4.28 fm. The potential giving the correct wavefunction 
has Vo = —19.54 MeV and a = 2.96 fm. Treating the Coulomb interaction as a first order perturbation about that 
finite spherical well interaction, wc estimate the Coulomb energy to be (e^/?-) = 0.57 MeV. 

The binding energies wc find for the nuclei with the chemical properties of helium (fUci), beryllium (2662), and 
carbon (3C3) by considering the analog nuclei and correcting for the Coulomb term are, respectively: 1.7 MeV, 23.9 
MeV, and 22.7 MeV. Thus, in the universe of light d and s quarks, the decay of carbon by fission into beryllium and 
helium is exothermic by 2.9 MeV. Since |Be2 has the same nuclear shell structure as the a particle in our universe 
(|He2), this fission reaction is the equivalent of a decay, which now occurs for nuclei as light as A = 6. The underlying 



Since d and s have the same charge, one might worry that weak transitions between the two participating baryons are flavor-changing 
neutral processes, and therefore suppressed by the Glashow-Iliopoulos-Maiani (GIM) mechanism. In our world, however, the heavier 
S~ decays into a S~ and a photon with a partial lifetime ~ 10^^* s [Tsll . primarily via a "long-distance penguin" diagram. This can 
be (very roughly) thought of as the tree-level process S~ — > A -|- tt~ , followed by the strong recombination of the it~ and A, to form a 
S~ (with a photon emitted). Here we will assume that, when energetically favorable, the conversion between S~ and E~ can always 
occur over time scales short compared to those needed for the evolution of intelligent life. 
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physical reason for this instabiUty can be understood in our world: the relatively weak nuclear binding energy of ^Lig 
compared with |He2 is almost enough to make it unstable to decay into deuterium and an a particle. The additional 
Coulomb energy in the (—1/3, —1/3) world is enough to generate instability. 

Through the SEMF and consideration of some specific analog nuclei, we find that heavier nuclei also undergo a 
decay, a process which is exothermic by about 8 MeV. In our universe a decay may be impeded for heavy nuclei by the 
Coulomb barrier through which the a particle must tunnel. In this case, however, the Gamow model of a decay (50l | 
shows that chemical carbon would have a very short lifetime — of the order of 10"-'^^ seconds. 

Universes in which the two light quarks both have charge —1/3 do not have stable carbon (or, for that matter, any 
nuclei with charge greater than four) and are therefore uncongenial. For universes with two light quarks of charge 
2/3, such as the u and c quarks, even more instabilities occur. Both participating baryons have a charge of 2, so 
the charge goes as 2 A. The two single baryons now correspond to isotopes of chemical helium (jHeg and oHC]^), and 
the stable a particle is chemical oxygen. The lack of hydrogen necessarily makes this an uncongenial world. As for 
carbon, there are two isotopes, ^Cj^ and f C2 (analog to fHe and fH in our world), but the additional Coulomb energy 
in this world destabilizes them. 



C. Two light quarks of different charges 

1. Calculations with analog nuclei 

Our universe has two quarks that arc light enough to participate in nuclear physics, u and d, with charge +2/3 and 
— 1/3 respectively. For all such worlds wc will revert back to the usual notation of N as the number of neutrons and 
Z as the number of protons in a nucleus, denoted by "^El. 

Let us consider first the case of Mp > Mn- Examination of the various isotopes of carbon in our world indicates 
that ^^C is the last to remain stable as Mp — M„ increases. The ^"'C nucleus in our world has a binding energy of 
105.3 MeV. It is stable against a emission (which is endothermic by 8.35 MeV) and other fission processes. It would 
undergo electron capture to form ^"^B if Mp — A/„ > 19.34 MeV, and decay by weak neutron emission (see Eq. (jIV.30|) ) 
into if Mp -Mn> 20.32 McV. 

The cutoff for congeniality, however, is not set by the requirement of stable carbon, but rather by stable hydrogen: 
When Mn < Mp + rrie, free protons may decay into neutrons by electron capture. Deuterium could serve as a hydrogen 
isotope, but it is only weakly bound (by 2.22 MeV) and easily undergoes weak neutron emission if Alp — M„ > 1.71 
MeV. 

In our universe, the /3-decay of tritium ('^H) into ^He is exothermic by about 0.02 MeV because, even though 
tritium is bound more strongly than ■^Hc by 0.76 MeV, the neutron is heavier than the proton by 1.29 MeV. With a 
proton-neutron mass difference large enough for free protons to decay, tritium would not undergo /3-decay. Instead, it 
could only decay through weak neutron emission and subsequent dissociation into three free neutrons and a neutrino, 

^H + e" ^ 3{n) + i^e ■ (V.5) 

The congeniality cutoff set by requiring some stable hydrogen isotope is therefore 

Mp - Mn < B{^R) - m, = 7.97 MeV . (V.6) 

This is the main result for the limit of congeniality for the case Mp > Af„.^^ 

Let us now consider the case Af„ > Alp. In this case, congeniality is determined by the stability of chemical carbon 
against all possible decay pathways. Whether or not a decay and fission would be exothermic depends only on the 
binding energies of the nuclei involved and not on quark mass differences. In our world, ^^C has a binding energy of 
92 MeV and is stable even against fission into three a particles. ^"C has a binding energy of 60.34 MeV, making it 
stable against fission into two a particles and two free protons by 3.7 MeV, and thus even more stable against simple 
a decay (which would leave an unstable ^Be behind). 

Weak proton emission changes a given carbon nucleus into a lighter isotope. ^^C decays by weak proton emission 
into "C ff AIn - Alp > 19.24 MeV; decays if M„ - Alp > 13.64 MeV; decays if M„ - Adp > 21.80 MeV, and 
decays if Af„ - Alp > 14.77 MeV. 



Even though wc have ehoscn not to deal with issues of nucleosynthesis in this paper, it is worth repeating the observation made in Sec. 
II Al that stellar burning might be very different in this world, since the two tritium nuclei can make ^He plus two free neutrons by a 
strong interaction. 
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Finally, /3-decay occurs spontaneously if 

M{Z,N) > M{Z + 1,N -1) +me . (V.7) 

The various carbon isotopes /3-decay into nitrogen if Af„ — Mp is greater than the following cutoffs: 18.64 MeV for 
i^C, 15.25 MeV for "C, and 24.40 MeV for i°C. Thus, ^^C is stable against all decays if M„ - Mp < 18.64 MeV. 

becomes stable against /3+-decay already if M„ — Mp is greater than a few MeV, and remains stable until it 
weak-proton-emits into for M„ — Mp > 21.80 MeV. For such large neutron-proton mass difference, the nucleus 
also weakly-emits into which then would /3-decay away. Thus, the congeniality cutoff for the case of Af„ > Mp is 
set by the stability of ^"C against weak proton emission: 

Mn - Mp < 21.80 MeV . (V.8) 

Even though we have not included the requirement of stable oxygen in our definition of congeniality, as mentioned 
in Sec. II Al most worlds with stable carbon and hydrogen contain also stable oxygen, and therefore the worlds we 
identify as congenial would usually allow the formation of water. This world is no exception. 



2. Estimation of dB /dA from the SEMF 



In analyzing the boundaries of congeniality for worlds with two light baryons, we have been able to use the analog 
nucleus method to determine the stability of carbon, without ever referring to the generalized SEMF constructed in 
Sec. IIV Bl However, it is instructive to study the SEMF predictions about stability boundaries against weak proton 
emission for different nuclei. These are shown in Fig. [Ml together with the predictions of the analog nucleus analysis. 
The bounds on M„ — Mp are qualitatively similar to the bounds we obtained from study of carbon and oxygen by 
the analog method. However, given the accuracy of the SEMF's predictions for femax(^), it is at first surprising that 
the SEMF's predictions of dB/dA w 5M/2 are off by as much as - 2 MeV. 

To understand why the SEMF docs poorly in this case, one must bear in mind that the parameters of the SEMF — as 
quoted in Table ITVl — are chosen so that the formula provides the best fit to the binding energy per nucleon B/A. We 
may write 

| = Ty(yl) + M^error(A) . (V.9) 

Note that the Weizsaeker formula provides an excellent fit for A > 10, with |Wcrror(^)| < 0.2 MeV. However, 

^ = W{A) + A^^^ + W^e,,or(A) -f (M/e„.or(A)) . (V.IO) 

Since the magnitude of the variation in Wcrror from A to A -f- 1 is genericallly of order of Wcrror(^), the last term 
in Eq. (jV.lOp can be quite large. If we were to use the SEMF to estimate the boundaries of stability against weak 
neucleon emission, it might be necessary to refit the parameters to optimize the estimate of dB/dA, rather than B/A. 



3. Congeniality bounds on u, d masses 

The condition for congeniality we found above for universes with light u, d quarks is 

- 7.97 MeV < M„ - Mp < 21.80 MeV . (V.ll) 

To translate this into a condition on quark masses, we could use the results from SU (3) perturbation theory presented 
in Sec. IIIII However, note that the best fit to the parameters of the model gave us an error of about 2 MeV in the 
masses of the proton and the neutron. For most of our purposes this is negligible, but in the case of worlds with 
only light u and d quarks we can improve on that accuracy by considering isospin symmetry. Taking into account 
electromagnetic corrections, the p-n mass difference is then 

Mp - Mn = (1.75 + 0.63 + 0.13) MeV . (V.12) 

Therefore Eq. (jV.lip translates to 



- 12.9 < xz< 4.1 



(V.13) 
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FIG. 14: Upper limits on M„ — Mp from stability against weak proton emission of nuclei at the bottom of the valley of stability, (a) for 
decay transitions involving a pairing term penalty, and (b) for decay transitions in which the nuclear binding energy gains a pairing bonus. 
The black curve is the prediction of the generalized SEMF, while the blue and red datapoints are obtained by considering analog nuclei. 



Referring back to Eq. (jIII.13[) . and taking m® = 100 MeV, this translates into limits, 

- 22.3 MeV< TO„ -TOd < 7.1 MeV . (V.14) 

Our universe with Mp — Mn = —1.29 MeV is at X3 — —1.17, comfortably away from the edges of the congeniality 
band. 



4- Nuclei with more than two baryon species 

If the s quark becomes light enough, a third baryon, the Ciow (which in our world is the A) will start to participate 
in nuclei. This marks the boundary along the xs direction of the congeniality region for two light u, d quarks: while 
universes beyond this boundary may still be congenial, the analysis in Sec. IV C 31 no longer applies. 

Consider the case X3 = 0, so that Af„ = Mp. Our universe is almost indistinguishable from this one in terms of 
nuclear composition. In the following analysis we will refer to the Ciow simply as A, both for the sake of the reader's 
intuition and because if X3 = 0, the Hamiltonian commutes with total isospin, so Ciow = A. 

Inside a nucleus composed of Z protons and N neutrons, a regular nucleon (which without loss of generality we 
will take to be a neutron) will spontaneously turn into a A particle if 



M(Z, AT, TVa = 0) > M(Z,Ar- 1, A^A = 1) 



(V.15) 



From now on when the mass M or the binding energy B is given as a function of three particle number variables, the 
last entry will be A^a (the number of A's). 

Hypernuclci have been experimentally investigated in our world [Tsj . For a hypernucleus that contains only one A 
particle, the binding energy of the A as a function of Z and A'^ can be expressed as: 



Sa(Z, N) = B{Z, AT, 1) - B{Z, N, 0) . 

Thus the mass of a hypernucleus is 

M{Z, A^ - 1, 1) = Ma + M{Z, A^ - 1, 0) - Ba{Z, N - 1) 



(V.16) 



(V.17) 



Figure 2 of [51[ provides experimental values of Ba{Z, N) for various hypernuclei. For large nuclei around A = 120, 
the A binding energy roughly saturates to _Ba = 23 MeV. For normal nuclei in our universe, B/A saturates to around 
8 MeV, so 

M{Z, AT - 1, 0) = M{Z, N, 0) - Mn + 8 MeV . (V.18) 
Combining this equation with Eq. ljV.lSP and Eq. (|V.17[) we find that the A participates in nuclear physics if 

Ma - Mn < 15 MeV . (V.19) 
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FIG. 15: Congeniality results on the triangle slice of xt = = 100. Our universe is the point by the bottom right eorner marked "us." 
The green bands are congenial, the red background is uncongenial and the central white region corresponds to worlds with more than two 
baryons participating in nuclear physics, worlds which we cannot at the moment fully characterize as congenial or uncongenial. Fuzzy 
borders imply uncertainty in determining boundaries. 

A theoretical calculation using only the kinetic Fermi gas model energy yields a cutoff of M\ — Mn < = 
24 MeV. The discrepancy between this result and the value of 15 MeV deduced from real hypernuclei must arise from 
A-nucleon interactions, which are not included in the two flavor SEMF. These interactions discourage the A from 
participating in nuclei until the A — N mass difference decreases to ~ 15 MeV. Also note that performing the previous 
analysis for carbon, we find that the A will not participate in carbon nuclei unless A/a — < 8 MeV. Therefore, we 
do expect universes to remain congenial a little bit beyond the boundary defined by Eq. (jV.19[) . 

This analysis applies rigorously only around X3 = 0. In that case the constraint on xs for the A not to participate 
in nuclear physics is 

X8 < -3.5 . (V.20) 

For l^sl > , the boundary between the worlds with nuclei containing only two baryons and worlds with nuclei 
containing three or more baryons would lie at xg < —3.5. The reason for this is that, due to the asymmetry energy 
penalty for nuclei composed of the two lightest baryons when these are not degenerate, it can be energetically favorable 
for a nucleus to capture a heavier A. 

Figure [15] shows the quark mass triangle for xt = 100, with the areas corresponding to congenial worlds colored 
in green and the areas corresponding to uncongenial worlds colored in red. Our world corresponds to the point in 
the lower right-hand corner labelled "us." The central region — where more than two baryons participate in nuclear 
physics — is shown in white, with a red dot in the middle to indicate that worlds where all three quark masses are 
identical are very probably uncongenial, as we shall argue in Sec. IV E 21 The rationale for the two upper, narrow 
green bands shall be discussed in Sec. IV Dl 

D. One light quark with two participating baryons 

If one quark is much lighter than the other two, then the two lightest octet baryons may participate in nuclei, 
leading to congenial domains different in character from the one in which wc live. These are the type (ii) regions 
shown in Fig. O They lie close to the midpoint of an edge of the quark mass triangle. To be definite, we consider the 
case where <^ rrig w niu. The other charge possibilities are considered later. For md small enough and rriu and 
close enough, there is reason to expect nuclei composed of neutrons and to form and to have qualitatively similar 
systematics to nuclei in our world. This region corresponds to P « 1/6 in the perimeter plots shown in Figs. [TUI and 

m 

On the other hand, there are important differences that distinguish this congenial region from type (iii) regions like 
the one in which our world lies. Type (ii) regions are smaller in extent, and harder to characterize quantitatively. 
First, as discussed in Sec. lV Al if to„ and are too large, the decuplet state (in this case A~ or ddd) may undercut 
all octet baryons, leading to an uncongenial world. If this happens for a given XT^ it excludes a region close to the 
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FIG. 16: Congeniality results on the triangle slice of (a) xj- = 50 and (b) = 200. The coloring conventions are the same as for Fig. 
1151 Both (a) and (b) arc drawn on the same 2:3, g scale, but not on the same scale as the triangle in Fig. 1151 



boundary of the quark mass triangle near P = 1/6. For example, in the case treated at the beginning of this section 
(to„ = 0) we found that the A++ became the lightest baryon if xt > 122 ± 29 MeV. Note that at the lower limit of 
the uncertainty this would occur at the perimeter of the quark mass triangle with xt = 100, in which our world sits 
(see Fig.[T5|). 

As the light quark mass increases (moving into the interior of the quark mass triangle), this problem diminishes, 
but another one appears. Notice (see Fig. [TT|) that for the same xt, other baryons are much closer in mass to the 
lightest nearly degenerate pair in this region than in type (iii) regions. In the case of a light d quark near P = 1/6, 
the proton, S~, and Ciow have masses less than 75 MeV heavier than the neutron and S~. In contrast in the type 
(iii) congenial domain near P = 0, the A is more than 170 MeV heavier than the p-n doublet. As a result, other 
baryons begin to participate in nuclear physics sooner (as and are increased) than in our world. We make a 
rough guess of how far from the center of the mass triangle this congeniality band ends, by supposing that, as in the 
case of the light u and d worlds, the third lightest baryon must be less that 15 MeV heavier than the two lightest 
baryons for it to participate in nuclear physics. Translating this into quark masses we obtain the boundaries shown 
in Fig. [151 

Any attempt make quantitative estimates of nuclear binding in type (ii) worlds is obstructed by the fact that the 
pseudoscalar meson spectrum is qualitatively different from our world. This can be seen in Fig.[T31 Near P = 1/6 the 
hadron spectrum has an SU(2) symmetry similar to isospin. However the lightest pseudoscalar meson is an isosinglet 
similar to the rj in our world. Not much heavier is a complex doublet similar to the kaons in our world. We suspect 
that the dominant contribution to the nuclear force still comes from exchange of a scalar isoscalar meson analogous 
to the /o(600), but we have no tools to analyze the problem quantitatively. 

It is possible to say something about the width of the congeniality band for type (ii) congenial regions. The strength 
of nuclear binding sets a limit on the n — S~ mass difference similar to the limit we derived from considerations of 
weak nucleon emission in the previous section. This translates into a limit on the s-u quark mass difference. However 
the neutron-S^ mass difference grows twice as fast with m„ — than the proton-neutron mass difference grows with 
niu — rud- All else being equal, therefore, the congeniality band has half the width of the one described by Eq. (jV.lSp 
in which our world sits. This is shown in Fig. 1151 

Finally we mention the extension of these ideas to other quark charges. First consider worlds where the three light 
quarks have charges 2/3, —1/3, and —1/3 (uds). We have discussed the case where the d quark is very light. A light 
s quark yields identical results. In worlds where the u quark is light, the proton (uud) and (uus) are the possible 
constituents of nuclei. However since both have positive charge, we do not expect congeniality, for the same reason 
that we failed to find a congenial region when the E~ and S~ were the light doublet. Finally, if the light quarks have 
charges 2/3,2/3 and —1/3 {cud), the candidate worlds with one light quark have either two neutral baryons, when 
the d quark is lightest, or baryons of charge +1 and +2, when the c or w quark is lightest. Neither can be congenial. 

For the uds worlds, the regions of congeniality are shown in Fig. [12] for the quark mass triangles with xt = 50 and 
with XT = 200. Notice that the width of the bands of congeniality, and the size of the central region where more than 
two baryons participate in nuclear physics, are independent of xt- For large xt, however, the narrow upper bands 
of congeniality (corresponding to neutron-E^ worlds) end before reaching the perimeter, for the reasons discussed in 
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Sec. |V3 

E. Three light quarks 

There are two distinct ways in which more than two nucleons could participate in nuclear physics. First, more than 
two quarks could have nearly degenerate masses. For three light quarks, this would correspond to area very near the 
center of the quark mass triangle, for any value of xt- However, like the accidental degeneracy discussed in Sec. IIIB| 
this would be quite unlikely if the a priori quark mass distributions are independent of one another. The second and 
more relevant case is that in which three (or more) quarks are light even though their mass differences might be of the 
same order as the masses themselves. Both cases can be treated simultaneously because the essential issue is flavor 
S'C/(3) breaking. 

If quark masses are logarithmically distributed (as we have argued they might well be) then the landscape would be 
dominated by worlds in which many of the quark masses are very light compared to Aqcd • If we add leptons to the 
discussion, we already have evidence from the neutrino sector of the possibility of very small masses. Unfortunately, 
the analysis of worlds with more than two baryons participating in nuclear physics requires better methods than those 
currently available. 

It is possible, nonetheless, to make some qualitative statements about worlds with more than two very light quarks. 
To simplify matters we consider the case where the quark masses are negligible and the only contributions to the 
baryon mass differences are electromagnetic. Since we are limited to considering at most three light quarks, there 
are four charge assignments of interest: (a) 2/3, —1/3, —1/3 (which we label uds for convenience); (b) 2/3, 2/3, —1/3 
{cud); (c) -1/3,-1/3,-1/3 (ds6); and (d) 2/3,2/3,2/3 [uct). 

1. A possible dihyperon instability 

The first relevant question about the domain of three light quarks is whether octet baryons are stable. It is possible 
that in this limit the dihyperon, H — composed of two quarks of each flavor — could be less massive than two octet 
baryons [1^1. If this were the case, then regular nuclei would decay into a state made of A/2 dihyperons. Since the H 
is a boson and therefore not subject to the Pauli exclusion principle, this would lead to worlds radically different from 
ours. Because no state made of dihyperons has unit charge, organic chemistry as we know it would not be possible. 

The various unsuccessful searches for a stable H suggest that, in our world, H is close to or heavier than twice the 
mass of A, 2,231 MeV, which is the threshold for its strong decay. For Mh < 2Ma, the dihyperon would still not 
be perfectly stable, since it could decay slowly by strangeness-changing weak interactions, unless it were also lighter 
than twice the nucleon mass (w 1,880 MeV). Thus 2Mn sets the absolute threshold for the stability of H . 

The quark content of the H (uuddss) matches that of two A's, two E's, or TVS. Flavor SU{3) violation lowers the 
mass of two As relative to EE or TVS. The H, on the other hand, is an SU (3) singlet and its mass is independent of 
SU{3) breaking. Since SU{3) becomes exact when all three quarks are massless, a measure of the stability of the H 
in that limit can be obtained by comparing its mass in our world with twice M (the average octet baryon mass) or 
2 X (1, 151) = 2, 302 MeV. Thus, in the absence of SU (3) violation, the threshold for H stability with respect to strong 
decay would rise from 2Ma = 2,231 MeV to 2AI — 2,302 MeV. The threshold for absolute stability is even more 
dramatically affected: It is ~ 1,880 MeV in our world, but would be 2,302 MeV without 5/7(3) symmetry violation. 

This argument is far from rigorous, but it does suggest that the H might be more stable in a world of massless 
quarks than it is in our world. This is an important point on which a lattice calculation might shed light. Determining 
the stability of _ff in a world of massless quarks would be interesting from an environmental perspective, especially 
if the qualitative features of nuclei built of dihyperons could be established. One might speculate, for example, that 
an instability to forming dihyperons puts an environmental lower bound on the number of very light quarks in a 
congenial world, in which case we might owe our existence to (among other things) the relatively large mass of the 
strange quark in our world. For the rest of this section, we assume that the dihyperon is not stable when three quarks 
are very light. 

2. Worlds with massless u, d, and s quarks 

Of the four charge assignments listed at the beginning of this section, the first, 2/3, —1/3, —1/3 is special because 
the sum of the charges is zero. It is also the one of most immediate interest to us, because it is closest to our world 
and occupies the apex of the quark mass prism we have been exploring. 
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If all octet baryon masses were exactly equal, then the minimum of the valley of stability would occur at electric 
charge neutrality, i.e., at equal numbers of each quark species. At low A, charged nuclei would quickly /3-decay or 
electron-capture until they reached the stable, electrically neutral nucleus with baryon number A. The only possibility 
of a stable, charged nucleus would develop at very large A, where the asymmetry term in the SEMF becomes less 
important. Thus, even small electromagnetic mass differences make a big difference in worlds where u, d, and s quarks 
are all massless. To estimate the effect, we used the electromagnetic mass shifts given in Table HID together with 
the three-flavor SEMF of Eq. pV.28p (in which baryon species abundances have already been adjusted to maximize 
binding at fixed Y and The parameter that drives the appearance of stable charged nuclei for zero quark masses 
is AM/ ~ —0.15 MeV. We then compute the lowest value of A at which carbon {Z = 6) is a stable against /?-decay 
and/or electron capture. First we minimize the nuclear mass with respect to Y since weak non-leptonic interactions 
will effectively adjust the d and s quark abundances to maximize binding. This yields ymin as a function of A, I^, 
and the SEMF parameters. Then we re-express the nuclear mass in terms oi Z = I-^ + ^Y. We focus on Z = 6. For 
small values of A, the state with Z = 6 decays to lower charge by electron capture. The criterion for it to be stable 
is therefore 

M{Z = 6,A)+me> M{Z = 5,A) . (V.21) 

The value of A at which this is satisfied depends most sensitively on the strength of the asymmetry term in the SEMF, 
Ca.^^ For €a = E® wc find that carbon becomes stable only for yl > 2 x 10^. Even if such nuclei were stable against 
particle emission, the likelihood of these gargantuan atoms coming together to form life forms seems very remote. 

The existence of very light, strongly bound nuclei analogous to the a particle in our world, would very likely prevent 
the formation of large A nuclei in worlds with massless u, d, and s quarks. Shell effects and the pairing term both 
favor the binding of nuclei made of two baryons of each species, starting at A ~ 4 (in direct analogy to the a particle) 
and continuing up to A = 16, a particle we call the "arkon^\ after its Biblical analog which carried every species 
two-by-two. These s-state nuclei would have their own complicated stability pattern, but there is every reason to 
expect one or more of these species to be very tightly bound relative to generic heavy nuclei (as is the a particle in 
our world), as well as electrically neutral. Heavy nuclei like the A « 10^ carbon which we have just considered might 
then very well be unstable against emitting such light, neutral objects. 

We conclude that there is good reason to believe that a world with very light u, d, and s quarks would be uncongenial. 
Quantitative bounds on the quark masses, however, await better theoretical tools. 

3. Other light quark charge assignments 

The three other charge assignments listed at the beginning of this section: cud, dsb, and uct, are intriguing. In the 
last case, uct, all octet baryons have charge Z = +2, so there is no hydrogen and these worlds are uncongenial. The 
minimum of the nuclear asymmetry energy occurs at Z = +A for cud worlds and Z = —A for dsb worlds. This implies 
a greater mismatch between the Coulomb and asymmetry energy minima than in our world (where the asymmetry 
energy minimum occurs at Z = A/2). Since the confiict between asymmetry and Coulomb energies destabilizes 
large A nuclei, we expect that, given similar SEMF parameters, the nuclear mass table in these worlds would likely 
terminate at smaller values of A than in our world. On the other hand, the very light nuclei that are important to 
us, hydrogen, carbon, and oxygen, might well be stable. For example, in the cud world the baryons with the quark 
content ecu and uuc have Z = +2, so a baryon number four nucleus with the structure of the a particle would have 
the chemistry of oxygen. Intriguing though these worlds may be, we leave them for future study, hopefully with more 
powerful analytic tools. 

F. Final congeniality results in a light u, d, s world 

We can now complete the description of the congeniality triangle introduced at the beginning of this section. The 
results are summarized in Figs. [15] and [16] . The green congeniality bands beginning in the corners at P = and 
P = 2/3 are the type (iii) congeniality domains discussed extensively in Sec. IV Cl The bands beginning at the centers 
of edges, at P = 1/6 and P = 1/2 are the type (ii) congeniality bands discussed in Section fVPI The center of the 



Weak scattering processes like An ^ nn are mediated by us <-> du at the quark level. 

It also depends on the Coulomb interaction, but we do not see any reason to change ec from its value in our world. 
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triangle is red because the neutral worlds discussed in this subsection are very likely to be uncongenial. The domains 
shown in white are regions where three or more baryons participate in nucleus building, but the quark mass differences 
are too large to apply the symmetry arguments of Sec. IVE2I Note that the type (lit) congeniality domains extend 
further into the triangle than the type (ii) domains because the next lightest baryon is split further from the lightest 
doublet in that case. 
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APPENDIX A: FIRST ORDER PERTURBATION THEORY FOR BARYON MASS DIFFERENCES 

Even though it is a textbook exercise, for the sake of completeness we present here the first order perturbation 
analysis relating octet baryon mass differences to quark mass differences. We label states in the baryon octet by 
|8; /a, Y, v) = |8; /i) , where v is the label that distinguished states of the same /s and Y that mix due to SU (3) flavor 
violation. The mass of an octet baryon identified by \x is (8; \x\ H |8; ^). Using the Wigner-Eckart Theorem, we may 
express the expectation value of the last two terms in -fffiavor in Eq. (jllLSp as: 



(8;A^|e« |8;a^) 



8 8 8(^) 



(8||e«||8). 



(A.l) 



The double-barred term is the Wigner reduced invariant matrix element (WRME) and the term in parenthesis is 
the SU{2i) Clebsch-Gordan coefficient. The sum over the index 7 reflects the fact that there are two irreducible 
8-dimensional representations in the decomposition of 8 (g) 8 = 27 © 10 © 10 © 81 82 © 1. The 5/7(3) Clebsch-Gordan 
coefficient can be decomposed as 



/i a /i 



(^) ^ 



, Y J Y J Y 

^fi-'-^ ^a^a ^^-^^ 



(A.2) 



where the C symbol is the SU{2) coefficient and the term in parenthesis on the right-hand side of Eq. (|A.2[) is the 
isoscalar coefficient. Their values are quoted, for instance, in |52| (whose notation we have adopted). 

We may therefore express the masses of the octet baryons in terms of five parameters: a term Ao = iJqcD +™o (©o) 
common to all the baryons, the two WRME's, 771,3, Siiid rrig. For example, the mass of the proton can be expressed as 



TO3 
a/3 

+ ms 



^(8||e«||8)^ + i(8||e«||8), 
-^(8||e«||8), + i(8||e«||8), 



(A.3) 



where the SU{2) Clebsch-Gordan coefficients arc shown in front of the terms in square brackets. 
Following the notation of [53|, we define two new parameters proportional to the WRME's, 



so that Eq. (jA.3P becomes 
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^^(8||e«||8), , F^^(8||e«||8), , 



Mp = Ao 



fZF -D 



7718 + {F + D)m3 . 



(A.4) 



(A.5) 



Analogous results for the other baryons on the periphery of the octet weight diagram are summarized in the first six 
entries in Table [D 
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The Hamiltonian of Eq. pil.ip is not diagonal in the (/a, Y^) basis. In other words, not all of the basis states in the 
octet representation are mass eigenstates. If we label the eigenstates of P the A (/ = 0) and the YP {I = 1), then in 
that basis (A, E°) and in the I3 ~ Y ^ subspace, 

n / ™ ™ \ 

+ ^^EM . (A.6) 




Defining A = y^m| + m|, the mass eigenstates in the (A, S") basis are, 



This mixing is negligible (for most purposes) in our world, since m^/ms — 0.02, and A = 1.0002 mg ~ mg. Therefore 
to a few percent, Ciow ~ A, and Chigh ~ However, this mixing is significant in worlds where 1713 and rris are 
comparable. In those cases total isospin no longer approximately commutes with the Hamiltonian, and the states 
Chigh or Clow cannot be identified with total isospin eigenstates. 
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